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Notation 



Cp - Pressure coefficient defined in eq. (3) 
a - Yaw angle relative to laboratory space 
0 - Pitch angle relative to laboratory space 

P t - Total pressure 
P $ - Static Pressure 
a rp - Yaw angle relative to probe axis 
4> - Pitch angle relative to probe axis 

Pp - Pressure indicated by the probe 
otp - Yaw setting of the probe 
<p p - Pitch setting of the probe 

Subscri pts 

I 

II relating to probe positions shown 



A relating to A type probe 

B relating to B type probe 

Lo lowest value in an array 
rp relative to the probe axis 
Up highest value in an array 



III ! in figure 3 



IV I 



1 



dotation of PENPTS 



DD - denominator in the subsequent line of the program 

HT - penetration height (Cp of a calibration surface) 

IR - a flag with a value of either 1 or 2 to identify first 

or second penetration point. This flag also controls 
return of values to calling program or subroutine 

IB - a flag indicating the location of the last rectangle 
in the band checked for a penetration point. It pre- 
vents repetition of the scanning in the 103 Do loop 

NX - number of X values in the calibration surface grid (rvalues of the 
calibration surface) 

NY - number of Y values in the calibration surface grid ( values of the 
calibration surface) 

X(I) - X values in the grid 

XM - value of X at point B (fig. 5) 

XR(I) - two values of X returned by the subroutine as results 

XS - currently calculated value of X which is the present 
result and is later stored in XR(I) 

Y(J) - Y values in the calibration grid 

YC - intermediate value used in the subsequent line of the program 

YG - the Y position of the penetrating line* 

Z(I.J)- Values of the Z = Cp calibration surface above grid points 

ZX - value of Z either at point A or point C (see figs. 5, 

6 and 7) 

ZM - value of Z at point B 



* It should be pointeci out here that the penetrating line 
is parallel to the X axis at a height HT = Cp above the 
X,Y plane and at a distance YG from the X axis on the X,Y plane 
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Notation of INTSCS 



AA 


“ 


constant defined in eq. 8 


AB 


- 


constant defined in eq. 9 


ARIN 


- 


ordinate at which scanning starts 


ARLO 


“ 


lowest ordinate of calibration surface (as a 
condition it must be identical for all calibration 
surfaces used) 


ARM1 


- 


ordinate measured relative to probe 1 (the probe 
belonging to the first curve Cp = const.) 


ARM2 




ordinate measured relative to probe 2 (the probe 
belonging to the second curve Cp = const) 


ARP1 


- 


setting ordinate of probe 1 


ARP2 


- 


setting ordinate of probe 2 


ARRES (1) 


- 


abscissa value returned to INTSCS from PENPTS 


ARUP 




highest ordinate of calibration surface (as a 
condition it must be identical for all calibration 
surfaces used) 


AR2R 


- 


ordinate relative to laboratory space 


AR2RJ 


- 


ordinate relative to laboratory space of previous 
scan 


BA 


- 


constant defined in eq. 8 


BB 


- 


constant defined in eq. 9 


DAR 


- 


ordinate step for scanning 


HT1 


- 


penetration height of calibration surface of probe 
1 (Cp of probe 1 ) 


HT2 


- 


penetration height of calibration surface of probe 
2 (Cp of probe 2) 


IEPS 




a flag which is equal to 2 when the calculation is 
carried out as shown in fig. 11, otherwise it is 
equal to 1. (This information is required in the 
main program and is returned to it) 
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ISL 

ISR 

K 

NT 

N2 

RES1(K) 

RES2(K) 

R1 

R2 

RAIL 

RAIR 

RAJL 

RAJR 

RBIL 

RBIR 

RBJL 

RBJR 

X1(I) 

X2(I) 



- a flag which is equal to 10 if one of the left 
hand side results returned from PENPTS is 1000.0 
(i.e., is not a penetration point). Otherwise 
it is equal to 1 

- a flag which is equal to 10 if one of the right 
hand side results returned from PENPTS is 1000.0 
(i.e., is not a penetration point). Otherwise 
it is equal to 1 

- a flag indicating whether the first or the second 
intersection is evaluated 

- number of abscissas in calibration surface grid 
(fig. 4) 



- number of Ordinates in calibration surface grid 
(fig- 4) 

- evaluated and returned abscissa of intersection 

- evaluated and returned ordinate of intersection 

- radius in fig. 11 for probe 1 

- radius in fig. 11 for probe 2 

- abscissa of point AIL in fig. 9 or 10 

- abscissa of point AIR in fig. 9 or 10 

- abscissa of point AJL in fig. 9 or 10 

- abscissa of point AJR in fig. 9 or 10 

- abscissa of point BIL in fig. 9 or 10 

- abscissa of point BIR in fig. 9 or 10 

- abscissa of point BJL in fig. 9 or 10 

- abscissa of point BJR in fig. 9 or 10 

- abscissas of calibration surface of probe 1 or 
its transformation as determined by calling 
statement in main program 

- abscissas of calibration surface of probe 2 or 
its transformation as determined by calling 
statement in main program 
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X3(I) - abscissas of calibration surface of 

transformation as determined by call 
in main program 

X4(I) - abscissas of calibration surface of 

transformation as determined by call 
in main program 



Y1(I) 


- ordinates 


match i ng 


xi(D 






Y2(I) 


- ordina 


tes 


matching 


X2(I ) 






Y3( I ) 


- ordina 


tes 


match i ng 


X3( I ) 






Y4( I) 


- ordina 


tes 


matchi ng 


X4(I) 






Z1(I,J) 


- values 


of 


Cp 

P rp 


for 


XI (I) 


and 


Yl( 


Z2(I,J) 


- values 


of 


Cp 

K rp 


for 


X2( I ) 


and 


Y2( 


Z3(I,J) 


- values 


of 


Cp 

K rp 


for 


X3( I ) 


and 


Y3( 


Z4 ( I , J ) 


- values 


of 


Cp 

P rp 


for 


X4(I) 


and 


Y4( 



probe 1 or its 
ing statement 



probe 2 or its 
ing statement 
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Notation of Main Program VDR 



A 

ALF 

ALFA 

ALFC 

ALFCN 

ALFD 

ALFDN 

AL FI 

ALF2 

ALF3 

ALF4 

AM 

CPA 

CPB 

CPC 

CPD 

EPSPS 

EPSPSG 

FALF 

FPHI 

IA50 

IA60 

I CO PS 

IEPS 



variables in library routine IXCLOK used to evaluate 
computation times 

values of a returned to VDP from INTSCS 
a input when experiment simulation is carried out 

“rpIII 

new value of ^ r pjjj f° r next iteration 
^rpIV 

new value of « r pjy for nex t iteration 

“i 

a II 

a I 

“III 

mach number 

Cp rpl 

Cp rpl I 
Cp rpl 1 1 
^ P rp I V 

relative difference between present static pressure 
and its value in the previous iteration 

convergence criterion (on the static pressure) 

final a value (at convergence) 

final value (at convergence) 

a flag in IXCLOK 

a flag in IXCLOK 

number of PS corrections carried out to ensure the Cp's 
are in calibration range 

a flag governing the convergence criterion returned 
from INTSCS 
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I EPS! 

IEPS2 

I IT 
IS CAN 

IXCLOK 

NOCOPS 

NOSIM 

NX 

NY 

PDYN 
PHI 
PH IB 
PH 1 1 
PHI1 
PH 1 2 
PH 13 
PH 14 
PPA 
PPB 



- a flag governing the convergence criterion for probes 
I and II 

- a flag governing the convergence criterion for probes 
I and III 

- number of iterations on p $ 

- number of static pressure scans from initial static 
pressure guess upwards 

- system subroutine for computation time evaluation 

- a flag which when equal to 1 causes static pressure 
corrections to be carried out such that Cp's are 
always inside calibration range, and when equal to 2 
cause these corrections to be skipped 

- a flag which when equal to 1 causes the program to 
simulate velocity measurement experiments and when 
equal to 2 causes the program to reduce measured data 

■ number of calibration surface matrix abscissa's, 
a values (must be identical for all calibration 
surface matrices used) 

• number of calibration surface matrix ordinates, 

<f> values (must be identical for all calibration surface 
matrices used) 

• dynamic pressure 

• value of 4> returned to VDR from INTSCS 

VpII 

• <j> input when experiment simulation is carried out 



Vpl 
Vpl I 
Vpl 
Vpl 1 1 

signal of probe in position I 
signal of probe in position II 
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PPC 



PPC 

PS 

PS IN - 

PSN 

PSC 

PST 

PT 

RELXPS - 

VI RTIM - 

XA 

XAX 

XB 

XBX 

XPC 

XPD 

XRIN - 
XRLU - 

XRUP - 

YA 

YAX 

YB 

YBX 



signal of probe in position III 
signal of probe in position IV 
static pressure 

initial guess of static pressure 

new value of static pressure for next iteration 

corrected static pressure 

memorized corrected static pressure 

total pressure measured 

relaxation factor for convergence on the static 
pressure 

virtual computation time 

abscissas of calibration surface matrix of probe A 

transformed abscissas of calibration surface matrix 
of probe A 

abscissas of calibration surface matrix of probe B 

transformed abscissas of calibration surface matrix 
of probe B 

yaw setting of position III 
yaw setting of position IV 
a value at which scanning starts 

lowest a value of calibration surface (as a condition 
it must be identical for probes A and B) 

highest a value of calibration surface (as a 
condition it must be identical for probes A and B) 

ordinates of calibration surface matrix of probe A 

transformed ordinates of calibration surface matrix 
of probe A 

ordinates of calibration surface matrix of probe B 

transformed ordinates of calibration surface matrix 
of probe B 
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YPA 

YPB 

YRIN 

YRLO 

YRUP 

ZA 

ZAX 

ZAMAX 

ZAMIN 

ZB 

ZBX 

ZB MAX 

ZBMIN 



- pitch setting of probe A 

- pitch setting of probe B 

- 4> value at which scanning starts 

- lowest 4> value of calibration surface (as a 
condition it must be identical for probes A and B) 

- highest value of calibration surface (as a 
condition it must be identical for probes A and B) 

- Cp rp ^ values of probe A 

- transformed Cp ^ values of probe A 

- maximum of ZA array 

- mi nimum of ZA array 

- Cp^pg values of probe B 

- transformed Cp g values of probe B 

- maximum of ZB array 

- minimum of ZB array 
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Introduction 



Experimental knowledge of the flow field generated 
by rotating turbo impellers is of prime importance in the 
research and development of turbomachi nery . It is essential 
for the refinement of design methods, for the development 
of new flow models which include secondary flow and tip 
clearance effects, and particularly for the verification 
of new computer codes developed to calculate the flow through 
rotating blade rows. 

In recent years laser velocimeter techniques have 
been applied successfully to measure the flow both inside 
and downstream of rotors. (Ref. 1 for example). It has 
become clear however, that the laser techniques are only 
reliable in the hands of experienced investigators. A 
window which remains clean is essential, and seeding is 
usually required. Laser techniques do not measure the 
pressure field and usually can only measure two components 
of the velocity unless the axis of the laser is tilted. 
Difficulty is also encountered when measuring close to 
walls. Hence there are reasons to consider alternative 
techniques, particularly if they are simpler to apply routinely 
in stationary turbomachinery passages. Furthermore, the 
achievement of redundancy in measuring the flow field 
behind the impeller is itself a worthwhile goal. The present 
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work deals with the application of a particular system 
of small high response pressures probes at the exit of 
an impel 1 er. 

Measurements behind an impeller, in the stationary 
bladeless gap, are simpler to make than measurements within 
the rotating passages. Transducer probes can be installed 
through the stationary machine casing and the data transmitted 
without resort to slip rings or rotary transformers. The 
sensor is not subjected to the centrifugal field or to 
the vibration of the rotor. However, the flow to be measured 
is then fluctuating at blade-passing frequency and any 
system of sensors must be calibrated for a wide range of 
possible Mach number, pitch angle, yaw angle and pressure 
variation - and yet must be capable of the necessary fre- 
quency response. 

In Ref. 2, a method was described for using two semi- 
conductor pressure probes together with the technique of 
synchronized sampling, to obtain the distribution of the 
velocity vector downstream of a rotor. The geometries 
of the two probes, designated Type A and Type B, and their 
installation in the compressor annulus are shown in Fig 
1. It was argued that, in principle, by rotating the probes 
in yaw about their tips and controlling the sampling of 
the data from each probe to be at the same position in 
the rotor frame, the system of two separate probes could 
be used to acquire data at a point in the per iodic flow 



TYPE A PROBE 



TYPE B PROBE 
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from the rotor, corresponding to data normally obtained 
from the multiple sensors of four- or five hole pneumatic 
probes when measuring velocity in steady flows. The tech- 
nique promised the use of probes having the simplest geo- 
metry and thus avoided the large size, expense and unreli- 
ability of multiple sensor probes which incorporate multipl 
semi-conductor transducers (Ref. 3). Because of the simple 
sensor tip geometry (that of a cylinder at incidence to 
the flow) the unsteady response was likely to be as good 
as could be expected of any single physical sensor. 

The two-probe technique is strictly applicable only 
to pe ri o di c f 1 ows . Howe ve r, data obtained on successive 
rotations of the rotor can be averaged to eliminate fluctua 
tions which are not periodic. This was shown to be effect- 
ive in tests reported in Ref. 2 in which a single Type 
A probe was used to establish the peripheral bl ade- to-bl ade 
distribution of flow yaw angle. 

In order to obtain velocity from the pressure measure- 
ments which can be obtained from the two probes, the steady 
response characteristics must first be established in cali- 
bration tests carried out in a known, controlled, uniform 
flow. Second, a method of applying the calibration to 
measurements made in an unknown flow must be devised. 

In the present method two different approaches have been 
followed. In Ref. 4, a technique is given for representing 



and applying the probe calibration analytically. When first 
applied the method gave surprisingly good accuracy (1 - 2 %) since 
the method required that the probes had characteristics 
which could be well represented analytically. This in 
turn required that the probe tips be geometrically precise, 
a feature which was not present in the first generation of 
probes . 

A second, more general approach is reported here, wherein 
the calibration of each probe is represented by a two-dimen- 
sional array of pressure coefficients. The application of 
the calibration given in this form, in an unknown flow required 
the development of special numerical procedures. The purpose 
of the present report is to document the analysis and the 
Fortran program developed to apply the method. 

In its present form, the method does not require that 
the calibration "Surfaces" be symmetrical about any axis or 
be expressed in analytical form, but does require that the 
pressure coefficient be independent of Mach number. The lat- 
ter restriction could undoubtedly be removed by introducing 
additional iterative steps. Further, in the present method 
only five measurements have to be taken to determine uniquely 
a velocity vector at a point. Throughout the report, the 
Fortran program notation has been used to describe the physics 
and equations involved in the solution. 
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Mathematical model 



Assume the A and B type probes of ref. [1] (see also fig. 1) 
to be immersed in a three-dimensional steady flow field.* The 
pressure response of each of these probes in given gas is func- 
tionally described by four variables as: 



Pp = P p (<* • 4> > P T > p s ) 



( 1 ) 



If a pressure coefficient is defined as 



Cp 



rp 



P 

P 



P 

P 



s 



s 



( 2 ) 



The calibration surface of each probe is given in the general 
case in form of a matrix of values of Cp, where 



Cp = Cp (a , d) ) (3) 

K rp K rp rp T rp ; v ' 

The pressure coefficient defined in this way has only a second 
order dependence on Mach and Reynolds numbers in the range of 



* For our purpose, using the synchronized sampling, the flow 
field behind the impeller is steady, although the probes re- 
quire a high speed response because of fluctuations. 



1 5 



0 <Ma< 0.7 in turbulent flows (Ref. 2), so that, to first 
order their influence is neglected in writing eq . (3). 

If the type A probe is rotated about its axis into three 
different positions (i = I, III, IV) readings are taken, and 
the type B probe is fixed in position II and a single reading 
is taken, the following four equations can be written. 

CP Ai ’ V-~p' S ' Cp Aj (a,<t,) 

p _ p 

Cp BII ° P* ^ P $ “ Cp BI] (oi,<f) 

it should be pointed out here, to avoid misunderstanding that a 
and <(> are defined in a coordinate system relative to the 
machine axis and not relative to the probe axis. In the set 
of four equations (4) there are four unknown quantities, namely: 
a, <j>, Py and P^ . These are the quantities to be evaluated using 
the measured data. Together with the stagnation temperature they 
define the flow field uniquely. The four equations, resulting 
from the four measurements, should be sufficient to determine 
the four unknown quantities. 

But the problem is complicated by three facts: 

1) The calibration surfaces are not generally known in 
analytical form. 

2) The calibration surfaces are double valued in a and $ 
i.e., for a given Cp and a there exist two <j> 
values, or for a given Cp and <f> there are two a 
values which satisfy eq. 3. 
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3) As a result of the measurement the value of P p rather than 

of C is determined. 

P 

Since the calibration surfaces are not given in a simple analytical 

form the solution has to be numerical. An iterative procedure is 

required because P and not C values are measured. First P and 

P P s 

P-j. have to be guessed to yield values, knowing the four measured 
values of P , and using eq. A. The guess is then iteratively 
corrected to converge on the solution. However, convergence of the 
iterative procedure is complicated because of the double valued 
nature of the calibration surfaces. 

This method of iteration shown in Fig. 2 was attempted initially 
for the evaluation of a measured point. In practice the initially 
suggested procedure converged in some cases and diverged in others, 
depending on the values, and signs of a and <J> . This was not 
surprising as convergence on two variables is not likely to be a 
simple matter. However, in order that the measurement technique 
be useful, convergence on the correct solution for a general set of 
measurements, is absolutely necessary. In practice, this can 
certainly be achieved if one of the two iteration variables is 
obtained by measurement. Since the static pressure is a difficult 
quantity to measure even in a steady flow field, only the stagna- 
tion pressure measurement need be considered. It is possible that 
the time-varying stagnation pressure could be measured with a suit- 
ably designed Kiel probe. Data would then be taken by synchronized 
sampling from the fixed Kiel probe, from the type A probe 
rotated into two positions, and from the 




Fig. 2 : 



Unstable iterative solution using measurement 
of A type and B type probes only. 
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fixed type B probe. The method of solution is then as shown 



in Fig. 3 . 
conditions. 



The method shown in Fig. 3 proved to converge under all 
It is described in detail in the following pages. 




Fig. 3: Stable iterative solution using measurements 

of A type and B type probes as well as 
measurements of a Kiel probe. 
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Evaluation of a and (method) 



The best way to understand the evaluation of yaw and pitch 
is to look at it from a topographical point of view. Each of 
the type A and type B probes has a unique calibration surface 
Cp rp = Cp rp ^ a rp ’ Vp ^ ’ where “ rp and <f> r p are measured relative 
to the axis of symmetry through the sensor at the probe tip. 

The calibration surface in this representation is invariant to 
yaw and pitch of the probe axis relative to the laboratory 
space. The same calibration surface if represented in the form 
Cp = Cp (a , <f>) can be written as 

Cp = C p [ ( a p + a rp ), U p + * rp )] (5) 

where a and <j> are the probe tip axis angular settings 
P P 

relative to the laboratory space. It is clear from eq . (5) 

that Cp can be derived from Cp rp by a constant translation: 

a , <f on the a , <t> plane. As each Cp can be viewed 
p p T r rp 

as a hill with its peak at a rp = 0 and 4 > r p = 0. The Cp 
surfaces are the same hills with their peaks translated to 
a = ap and = <f>p . 

In the present method probe A is used in three different 
angular settings namely: 

!) a p = 0 f p = 0 
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Ill) 

IV) 



0 



a p a plll ^p 
a p = a pIV *p = 0 

and the B probe is used in a fixed position, namely: 



II) « p = 0 „ p 



♦pi I 



For this case the topography of the calibration surfaces will 
appear as four hills. The three with their peaks at points 
(0,0); C+ctpjjj, 0) and (-ct pI y, 0) are the translated hills Cp rpA 
and the fourth with its peak at ( 0 , 4 > p jj) Is the translated 
hill Cp rp g. Their contours of constant Cp are then as shown 
in figure 4. 

Assume now that a velocity vector with yaw a and pitch 
<j> is to be measured. These values of a and 4 will be 

sensed uniquely by the probes at their four angular settings. Were 
equation (4) single valued, the values of a and 4 could be 
uniquely evaluated, as the single intersection point between the 
projections of appropriate lines of constant Cp on the a , 4 plane*. 
In the present case of the double valued functions, the lines of 
constant Cp are closed curves and more than a single intersection 
point do exist. 



* Cp has generally a different value for each probe in each 
of its settings. Thus the solution involves solving for the 
intersection points between projections of contours of specified 
(but different) Cp values on the different hills. 
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Fig. 4: 



Projections of = const, lines of the four 

calibration surfaces (the center hill is for 
probe A at a p = 0 , 

for probe B at a 



H P 

= 0 , 



0, the top hill is 
*P 



= f j j , the ri ght 



hill is for probe A at ap = +ot p 1 1 1 9 4p = 0 

and the left hill is for probe A at a = a 

* = 0. P 



pi V ’ 
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The situation is shown in Fig. 5 for an example of a number of 
such intersection points. The correct intersection point, however, 
is uniquely identified as the only point through which all four 
Cp = const, curves pass. 

From an examination of Fig. 5 it is clear that the inter- 
section points of three closed Cp curves projections only are 
sufficient to identify a and <f> uniquely. However one of 
these three curves must be that belonging to the type B probe. 

The details of the numerical procedure used to obtain the 
correct intersection point in the evaluation of a. and $ , 

are given in the following paragraph. 
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Fig. 5: An example showing a number of intersection 

points. The real one is in the. second quad- 
rant. 
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Evaluation of the intersection point coordinates (numerical 
Procedure ) 

The calibration surfaces Cp rp are represented in form 
of a linear string of values ordered into an array, as shown 
schematically in figure 6. The numbers in Fiq. 6 indicate the 
position of a Cp rp value in the string. The strinq starts 
with the value of Cp rp at a point (-a rpLo ; -4 rpLo ) 
and ends with the value of Cp rp at a point ( +ot rpUp '’ + <J> rpUp )- 
This sequence must be kept and can not be changed. 

Each of the closed Cp = const, curves projections in figs. 

2 and 3 is the projection of the line of intersection between a 
plane parallel to the a , $ plane at a height equal to Cp and 
the calibration surface. These closed curves can be determined as 
the locus of the projections of all the penetration points of ar- 
bitrary lines parallel to the a plane at a height Cp and the cali- 
bration surface. Such penetration points are calculated in sub- 
routine "PENPTS", Fig. 10. 

PENPTS calculates the first two penetration points of a 
surface Z (in the present case Z = Cp) by a straight line 
piercing that surface. If the surface is double valued these 
two points are the only roots. The surface is given as a table 
of numbers on a cartesian basis Z = Z(X,Y) (or in the present 
case Cp = Cp( a , <(>) . 

The subroutine has the following limitations: 

1) No roots can be found on the lower Y = const, boundary 
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Fig. 6: The order of a calibration surface 

array . 
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2 ) 



3) 



4) 

5) 

6 ) 
7) 



The roots can be found only alone a piercing line 
which is parallel to the X axis (but at any height 
above the X ,Y plane) 

The table Z = Z(X,Y) defining the surface must be 
based on a grid comprising lines X = const, and 
Y = const. The spacing of these lines must not be 
equal. In other words the surface is defined by a 
rectangular grid in the X,Y plane, from ^ m -j n to 



X m3v ancl from Y • „ to 
max mi n 



max 



Only a single root can be evaluated in a surface 
element located above a defining rectangle on the 
X , Y plane. 

Not more than the first two roots will be evaluated for 
any pi e rc i ng line. 

The surface must be monotonic over each rectangle 
(this is a result of limitation 4). 

All X and Y arguments must be given in increasing order. 



These limitations do not restrict the application of PENPTS 
in the present problem as long as the calibration surfaces are 
smooth within the element located above a grid rectangle. 

However, the elemental grid rectangles can be reduced arbi- 
trarily in size. If a calibration surface is more than double 
valued PENPTS will fail. However in this case the probe yielding 
such a calibration surface can not be considered a useful instrument 
unless used only in parts of the domain where it is double valued. 
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Z is defined from X . to X , and from Y . to Y 

min max min max 

PENPTS is given the following initial information: X values and 

Y values defining grid points, corresponding Z values, the Y 
location (YG) and the height above the X,Y plane (HT) of the 
piercing line. With this information PENPTS searches for the band 
of rectangles which includes YG or of which YG is the lower bound- 
ary (see fig. 7) and then scans this band from left to right in 
search for penetration points. The scanning is based on the geo- 
metry given in fig. 7 which represents a particular element in 
the band, approximated by two plane triangles. 

Initially a coarse scan is carried out just to detect, but 
not to evaluate, an intersection point. This is done checking 
for each sub-domain whether (ZX(I-l) - HT)/(ZX(I) - HT) 0. 

If this condition is true a root is detected and control is 
tranferred to its exact evaluation. The value of the root is 
calculated after its location, either in the first (left) or 
second (right) triangle is determined (each grid rectangle is com- 
posed of two triangles). Equation 6 which is based on fig. 8 
(for a left triangle) or eq. 7 which is based on fig. 9 (for a 
right triangle) is used to evaluate the penetration point. These 
equations express the linear interpolation of Cp in the Fortran 
rotation used in this program. 

XS = X(I-l) + ABS ( Z X ( 1-1 )-HT)*(XM-X( 1-1 ) ) / ( AB S ( HT- ZM ) 4 
+ ABS ( ZX( I -1 ) - HT ) ) (6) 
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Y 




Fig. 7: The geometry of a linearized calibration 

surface element comprising two plane tri- 
angles and its intersection with a plane 
normal to X,Y along Y = YG. 
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Fig. 8: The geometry for a penetration through a left 

hand ( fi rs t) triangle. 
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Fig. 9: Tne geometry for a penetration through a 

right hand (second) triangle. 
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XS = XM + (X( I)-XM)*ABS(ZM-HT)/(ABS (ZM-HT) + 

ABS(HT-Z( I ) ) ) (7) 

Equations 6 and 7 are invariant to the slope of the cali- 
bration surface i . e . , the slope of the straight lines AB and 
BC. 

Finally it should be pointed out that when two penetration 
points are determined the values of their abscissas, X(I), are 
returned. When only a single penetration point is detected the 
abscissa X ( 2 ) will be returned with a value of 1000.0. When 
no penetration point is determined both X(I) values returned will 
have the value of 1000.0. The pro gram logic is designed to 
recognize these messages. 

It was stated earlier that PENPTS is used to determine the 
closed intersection curves projections on the X,Y plane. In fact 
not the curves but just the intersection points between each two 
of them are required (see figure 3). 

To compute the coordinates of these points the subroutine 
"INTSCS" is used. It uses PENPTS as a subroutine. 

In INTSCS a scanning procedure is carried out from a 
minimal value of Y (or 4 1 ) YRIN to an upper value of YRUP, or a 
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prescribed 10,000 times* which ever comes earlier. The sub- 
routine scans through any two arbitrarily chosen closed curve 
projections to find their intersection points. In each scan 
(I) up to four penetration points can be determined, while the 
penetration points of the previous scan (J) are memorized. 

Together, eight penetration points can be involved. When no 
intersection point exists the geometrical situation is as shown 
in figure 11, while the existance of an intersection point is 
characterized in figure 12. Subroutine INTSCS can distinguish 
between the two situations. In figures 11 and 12 the case of 
four penetration points found in each scan are shown. The 
subroutine, however will handle any possible number of such points, 
from zero to four. A "no penetration points" is assigned an 
abscissa value of 1000.0 by PENPTS as explained earlier. 

When an intersection point is detected its evaluation is 
based on the geometry shown in fig. 12. The eight penetration 
points have the following coordinates: 



* The 10,000 scans are governed by a program constant in the 
line "DO 140 ...." and can be arbitrarily varied. 



AJL( RAJL ,YRJ) 
BJR( RBJR.YRJ) 
B I L ( RB I L , Y R ) 

A I R ( RA I R , Y R) 



BJL( RBJL ,YRJ) 
A J R ( RA J R , Y R J ) 
AI L ( RA I L , Y R ) 

B I R ( RB I R , Y R ) 
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Fig. 11: 



The geometry for two successive scans 
when no intersection point exists. 
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Y=* 




Fig. 12: 



The geometry for two successive scans 
when an intersection point does exist. 
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The left intersection point illustrated in fin. 12 is the 
intersection point of the lines aa and bb, each described by 



its equation: 

for aa Y = ( AA ) X + BA (8) 

for bb Y = ( AB ) X + BB (9) 

The constants AA, BA, AB and BB are given in eouations 10 to 13. 
AA = ( Y RJ - YR ) / ( RA JL - RAIL) (10) 

BA = YR - AA*RA I L (11) 

AB - (YRJ-YR)/( RBJL-RBIL) (12) 

BB - YR - AB *RB 1 L (13) 

The coordinates of the intersection point to be calculated are 

X = (BB-BA)/ ( AA-AB ) (14) 

Y = X*AA + BA (15) 



These relations are true for a left hand intersection point. 
Analogous equations are true for a right hand intersection 
point. In this algorithm the straight lines aa and bb approxi- 
mate the curved lines connecting AJL with AIL and BJL with BIL 
or similar lines on the right hand side of figure 12. The error 
introduced through this approximation is reduced as DAR = A <j> 
is reduced. 

It is possible that an intersection point is identical with 
AIL and BIL or AIR and BIR. This case is defined as "direct hit". 
The program is designed to detect such a direct hit and evaluate 
the correspondi ng intersection point accordingly. 

The above algorithm works perfectly as long as the two 
closed Cp = const, curves are far from being tangent. But in 
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practice a situation of almost tangent curves can arise when X 
(or a) is very small. In this case the preceeding algorithm 
will fail and must be replaced. The geometry of this situa- 
tion is described in figure 13. This situation is identified 
by INTSCS and the intersection points are then evaluated assuming 
that they are intersections of two circular arcs. When the 
curves Cp = const, are almost circular this approximation does 
not lead to unacceptable errors. 

It was stated earlier that INTSCS scans from a minimum 
value of <f> to a maximum value of $ with prescribed steps 
A <)> , as shown in figs. 11 and 12. This direction of the scan- 
ning is used when the intersection points of the curves Cpj 
and Cpjj of fig. 5 are evaluated. 

However, in course of the reduction of the measured data, 
scans in the direction of a in steps of Aa are also neces- 
sary to evaluate the intersections of the curves Cpj and Cp^^ 
or Cpj and Cpjy. INTSCS is designed to carry out this task 
as well. To do this the calling statement for INTSCS is appro- 
priately changed as will be explained in the next section. To 
be general enough INTSCS is not written in terms of a and i or 
X and Y but rather in terms of general arguments. The best way 
to understand INTSCS is to compare its general arguments to physi- 
cal quantities by means of the calling statements. In fig. 14 the 
flow diagram in INTSCS is given in terms of the general arguments. 

INTSCS returns to the main program (fig. 15) the (a,<t>) 
coordinates of two intersection points between the Cp = const, 
curves specified in the calling statement. Let us now follow 
the way in which the main program is designed to utilize INTSCS 
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Y= 4> 




Fig. 13: The geometry when two closed C 

curves are almost tangent. 



■►x = a 



= const. 
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FIGURE 14. FLOW DIAGRAM CF INTbCS 



for the evaluation of a and 41 of the velocity vector, as 
well as static and total pressures. 
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Evaluation of the Velocity Vector and Pressures from the Probe 
Signals 

Fortran program VDR was written to evaluate the velocity 
vector from probe measurements of pressures. The program is 
shown in Fig. 15. 

At line 1410 INTSCS is called to scan curves I and II 
for possible intersection points. Scanning can be carried out 
in the direction of the ordinate only, with the calibration 
curve matrices compiled exactly as shown in fig. 6. This limita- 
tion is imposed by the way PENPTS is constructed. In the case 
of the intersection points between I and II ZA and ZB are 
scanned without difficulties in the direction of the ordinate 
which is <f> as shown in fig. 12 and returns with the coor- 
dinates of the two first intersection points, points 1 and 2 
of Fig. 16. They are ALF1, PH II and ALF2, PHI2. In line 
1500 INTSCS is called again to scan curves I and III. Now 
scanning has to be carried out in the direction of the abscissa, 
a task for which INTSCS was not designed. To overcome this 
problem the calibration curve matrices are used in a transformed 
form such that the previous abscissas are now ordinates, or- 
dinates are abscissas and the internal structures of the Cp rp 
arrays are accordingly modified. This transformation is carried 
in VDR in the section between lines 320 and 450. Comparison 
of the calling line 1500 to the previous calling line 
1410 shows very clearly how the various arrays: original and trans- 
formed, are used. The coordinates of points 3 and 4 of fig. 16 
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Fig. 16: Selection of the proper intersection 

point in VDR. 
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are now returned to the main proqram of VDR, they are ALF3,PHI3 
and AL F4 , PH I 4 . 

We are seeking a single intersection point; the one repre- 
senting yaw and pitch as sensed by probe A in positions I and III 
and probe B in position II. But, as the calibration surfaces 
are double valued we are now, unfortunately, in the possession of 
four points. The solution, however, is physically unique. Only 
a single velocity vector exists in reality and its yaw and pitch 
are included in the four intersection points evaluated. Were 
both the measurements and the numerical procedure absolutely 
accurate, two of the four points would have been identical. But 
this is not the case in reality, instead of a sincle point, two 
points close to each other will be detected. Therefore the 
average coordinates of the two of the four intersection points 
returned to VDR from INTSCS which are closest to each other 
are selected as the measured yaw and pitch angles. In the ex- 
ample shown in fig. 16 this will be the point dividing the dis- 
tance between points 2 and 3. 

The calculation is now at the point at which a and q. are 
temporarily known (see fig. 3). Using a and <J> a new 
C P r p I V P r °be A in position IV) is computed by linear interpol^ 

tion using subroutine INTPLT. With this new value of CP r pjv and 
with PPj V a new P $ is computed. If this new value of P s is close 
enough as determined by EPSPSG to the guessed value, or to the 
value of P $ in the previous iteration, the data re duct ion for the 
particular point in question is terminated. 
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The relative difference between present and previous P r 

s 

is compared to the converqence criterion EPSPSG. This criterion 
is evaluated by an empirical function determined to give best 
compromise between accuracy of results and computer time required 
until convergence is achieved. When the calculation follows 
normal routine the convergence criterion is given in line 
1710 as function of Py . When the routine for almost tangent 
curves is used during data reduction either on the right or the 
left side a different empirical function (lines 1720 or 1730) 
is used. 

If the convergence criterion is satisfied results are 
printed out and data for a new measurement point is read for 
reduction. If, however convergence is not reached a new static 
pressure for a next iteration is evaluated (line 1780). 
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Convergence and accuracy 



The convergence of the present iterative procedure is not 
ensured in all possible cases of data sets. The program, however, 
is adjusted to converge in most of the cases. Similarly to other 
iterative computation method some experience is required to achieve 
convergence when the calculation does not converge. In this para- 
graph the principal factors affecting convergence, computation 
time and accuracy, which are abviously coupled, are pointed out. 

The convergence and accuracy of reduction of a data set: Py, 

PPA, PPB, PPC, PPD, PS I N (or analogous set in the experiment 
simulation mode) depends on the following factors, which can be 
varied by the user. 



1) Coarseness of calibration arrays ZA, ZB and their 
linearity. With coarser arrays convergence problem 
will increase and accuracy of results decrease. 

2) Probe settings YPB, XPC, XPD (with the following probe 
settings being fixed and not variable YPA = YPC = 

YPD = 0.0; XPA = XPB = 0.0). The optimum setting is about 
^20° to 25°. Too small values reduce accuracy, too big 
values cause convergence problems and probe tip flow 
separation. 

3) RELXPS, the static pressure relaxation factor. The 
smaller this factor the safer will convergence be 
achieved. Computation time, however, will increase. 

4) The constant 5.0 in line 1050 of the program. The 
dimension of this constant is kg/m. If too big 
results can be lost and if too small computer time 
will be growing. The variation of this constant should 
be coupled with an appropriate modification of the 
constant in line 1010 of VDR (item 8 in this list). 
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5) The constants in the evaluation of EPSPSG (lines 1710 
to 1730). The smaller EPSPSG the more accuarte the 
results in expense of increased computation time and 
reduced convergence safety. 

6) The constant scanning step DAR. The smaller DAR the 
more accurate the results, but too small values can 
cause complete loss of results. Computation time in- 
creases with reduced DAR. 

7) The constant 10 in line 1360. This constant governs 
the number of P $ corrections. 

8) The constant 1000 in line 1010. This constant governs 
the number of scans. 



If convergence is not achieved in a particular case varia- 
tion of each or of a combination of the above values will always 
enable convergence. 
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Evaluation of the Computation Time 



By deleting the letters CT from column 1 and 2 of lines 
20 to 70, 1930 to 1960, 1980, 2000, 2010 the actual computation 
time as well as CPU time are evaluated and printed out. The 
following statement prior to execution is required in this case:* 
GLOBAL T SYSLIB SSPLIB 

This option is useful for adjusting the constants affecting 
convergence as to optimal compromise between accuracy, ease of 
convergence and calculation costs. 



* When the program is run on Naval Postgraduate School IBM 360 
system. 
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Input 



1) Calibration arrays have to be input in the following manner: 

first values of NX and NNY , second values of X in rising 
order, third values of Y in rising order and fourth values 
of Z in the order shown in Fig. 6. Two calibration arrays 
are read: first the array of type A probe and second the 

array of type B probe. 

2) If the program is run in experiment simulation mode PT 
and PS are read in Kg/m? and subsequently o and . 

3) If the program is run in data reduction mode the measured 
values of PT , PPA, PPB, PPC and PPD are read and then 
the guess of static pressure PSIN. All are read in Kg/m?. 

PS I N has to be lower than the actually existing value to 
ensure convergence of the calculation. Too low a value 
will cause waste of computer time. 

Calibration arrays are read from a disk space on which they 
are stored. The following statement prior to execution is re- 
quired for successful reading: 

FILEDEF 02 DSK NAME XX 

here NAME is the name of the file on which the calibration data 
of probes A and B is stored in proper order and format, XX 
is a two digit number. 

Simulation or measurements reduction input is read in the 
normal way using the terminal keyboard or punched cards. 



51 



Output 



a, <t>, P $ , Py and Ma are printed out according to lines 
1820 to 1920. This output, however is not sufficient when the 
logic of the computation is to be followed either to examine the 
execution of the program or for debugging. Two programs with 
additional output are given in appendices 2 and 3 and can be used 
for this purpose. The first, SWVDR gives short additional 
output, namely: 

1) When PSIN is successively increased automatically 
by the program to ensure convergence, the values of 
PSIN are pri nted out . 

2) 1 1 T , FALF, F P H I and PSN are printed out at the 
end of each iteration prior to the convergence test. 

The second program WVDR prints more detailed information. All 

additional WRITE statements in this program are numbered with 

three digit numbers starting with 9. 
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Cone! us ions 



The following conclusions are drawn on the basis of experi- 
ence gained in running the program with various data sets and 
various types of calibration surfaces. 



1) Optimal probe settings for ease of convergence and 
high accuracy are: YPB = 25°, XPC = 25°, XPD = -25°. 

It is therefore suggested that a B type probe with 
250 pitch will be used, and that the A type probe 

be rotated to +_ 25° . 

2) Convergence and accuracy, as well as computer time 
efficiency are improved when the calibration surfaces 
of both probes are not flat at their peaks but are 
rather rounded. It is therefore suggested that a new 
probe tip geometry be considered. A spherical tip with 
a central pressure tap is recommended. To prevent 
damage to the sensitive transducer located behind the 
pressure tap, and in order to improve the frequency 
response of the probes it is suggested that the volume 
between the pressure tap face and the transducer 

be filled with an appropriate liquid that will not 
affect the transducer negatively. In this case the 
opening of the pressure tap has to be sealed with a 
very thin low inertia membrane. 

3) It is probably possible to modify the iterative pro- 
cedure such that safe convergence can be achieved also 
when using the scheme of Fig. 2. If this can be 
achieved the Kiel probe will not be necessary. An 
effort in this direction is suggested. 
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APPENDIX I 



LISTING OF VDR 
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I o « <i 3 • 3-, 



or/oa/Y'j 



FILE: VDR 
C START IF VJR 



FORTRAN PI 



CT 


INTECf "i 


* (b J 


CT 


IAS 0=0 






CT 


I A c/l= L 








Ci^'CNS 1 


•VJ 


F ( iO 




1 L A 1 9 C , V 
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a * i<+ ) J 
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<- 1 -> J 
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1 J 
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F C h F lf-0. 
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F Ir'MATlFll 


• 


135 


FCH' 4 £ 7 ( F l'j 


.5 ) 


13c 


FJRYATl 


r 1 / 


• 5 i 



Y 3X t <t U) ,Z3X < 90 , -,0) 



C FEAO PRO ft CALI -i-.U INN D'->T\ 

Rt A J ( ? , l l 9 ) kA,NY 
00 iii 1=1, OX 

121 RtACfi, ,J)J \ 41 1 ) 

CJ 122 1 = 1 ,NY 

122 RE40(2, 1 3 C 1 rA( 1 J 

CO 1 23 1=1, l< 

yj 123 J = 1 , 1 Y 

123 REAC1 2, 1 2 J 1 2 U 1 ,J i 
Cl) 22 i 1 = i , V< 

22 1 REaD( 2, 1 30 A “ ( i 1 
CQ 222 l = i , IY 

222 FEAO 12 , 1 3C) Ylll) 

CO 223 1=1, 'IX 

CO 2 23 J = i , NY 

223 READ! 2 , 120) L j ( 1 , J) 

C TRANSFER M - AL I FRariJfi •'•A TRICES 
OL 229 1=1 « N 7 
229 XAX1 11= iA( 1 ) 

DO 22 1) 1=1, NX 
225 Y«x( i) =X4i i) 

00 2 26 1=1, 'll' 

EC 226 J=i , N< 

22 6 L AM 1 , J) =ZA( j ,1 ) 

00 227 1= 1 , 'IY 
22 7 XriX ( 1) =YP ( i ) 

00 226 1 = 1, 1 A 

228 Y E A ( I ) = A f ( 1) 

OC 229 1 =1 , \< 

00 2 2 9 J = 1 , 'I A 

229 Z‘3X ( 1, J )-LnU , 1 ) 

C FRCEt. ScTYlNcS 

YPA= C. ^ 

YP E= 55 . 

>P A= 0. J 
XPC= 20 . 

XP L = -2 J . 

C PROGRAF C L A S I A ,Y I 3 
YRL0 = -2J . 

YK L F =3 J . 

YR I N = - 3 C • 

XRLC=-jJ . 

>R LP =30 . 

>R IN=-jC. 

1CCFS=1 
NO CO P ? = I 
NOS I f, = i 
II T=1 

1 SC A N = 1 

R EL X PS= J .3 

C EXPE-il t>c 7 1 SI •'LL AT ION 
<,0 1 FOR^TIFI;.*,) 

193 FlAC 15, ,0 1) ->T 
luRlTEl 6 ,901 I PT 
IF (NCS 1 -1 . F j . ,. ) j I 70 500 
Rt AO l 5,401) Pi 



IN TOCO lu 
I .\ TO 0020 
IN T COC 30 
I J T0OJ9J 
INTO TO 90 
IITC JO 60 
INTO J 0 70 
INTCOCeO 
I NT 060 90 
IN 7 00 100 
IN T C C 1 1 0 
IN TOu l 20 
INTO ) 1 j ) 
IN TOC 190 
I N TOul 50 
INTCUcG 
INT0J1 7 0 
INTJJluu 
INTlC 1 SO 
INT0J2 J) 
I N T 0 U 2 1 0 
IN I CO 2 20 
I NT 0 02 30 
1 N T C C 2 9 0 
INTOOOiO 
I N T L 0260 
1 J Y 0 J 2 70 
I NT 002 <30 
INTCC2S0 
IN T CO 3 00 
1 NT 00 3 10 
INTCC320 
INI 003 30 
IN TC0340 
HTC0350 
INT0J36J 
I N T C 03 70 
I NT 0 0 380 
I NT 003 90 
MTC09C0 
1 NT Jj9 10 
I N T 069 2 J 
IN T CCA 2 0 
IN T 009 9 J 
INT0O990 
IN TC0960 
INT 0JS7j 
INTO^hBO 
I NT 00 990 
I NT JOS 00 
IN TC05 10 
IN C 0 5 2 J 
INTJ053J 
IN TCC59C 
INT 003 50 
I N T 0 w 5 o 0 
INT 00 6 70 
I N T 0 0 3 6 J 
INTCC550 
INTG0600 
1 N T 6 06 LO 
IN TCC6 20 
IN TCCfc 3 0 
INT0069J 
INTCU5C 
INTO Jt . 60 
IN TOOL 70 
IN T 0 n. 6 e 6 
INTO J690 
IN TCO 7 CO 



FILE: VJF FnPTPAM PI N 

READ! 5, 4 Ji) :LM “ 


0 s 


m jj 7 u 

I MOvi/^J 


l»Rl T ElA,tJl) A L r A 


, 


IN T C 0 1 i 0 


REALM 5, ^ C l J °Hi 1 




1N7 0 j74J 


wRITfc(o»tji) w hi 1 




IN T jJ / V) 


C4LL I (MlIUA ,Y a ,/A ,A LFA , PHI 1 , \X , “IY , CPA) 




I N i 0 J 7 o J 


alfC = a ..f\-ap; 




I NT 0 0/ / J 


CALL 1HT3L T (X A ,y A ,Z A , ALFC , PH I I , NX, ‘JY, C^C) 




1 N TC070u 


4 lF J-A lF A- a PJ 




IN f Co 7 J J 


CALL J 1 T ?L T (xA.y A ,Z A, A LTD, PHI I ,NX, iY ,CPJ) 




I N T j jd j.J 


PH 1 B= Fit I 1 -Y P d 




1 J TCCfi i J 


CALL 1 1 1 PLT U3.YP ,Zc , ALFA , PH I e ,NX , NY ,CPB) 




IN TOJ H<. J 


FD YN= PT -P j 




I *i 1 v# 0 i i J 


PPA -C P * rJ r\+PS 




1 I T C o H 4 J 


PP-3= CF’ l*PJYw +P 5 




in r j j j ^ j 


P FC=C PC * F j Y \ ■* PS 




iNTCCfcoO 


PPO = CP J> PJY Mt-Pb 




IN TO H 7 J 


C READ Me AStJPfc N LM T 3 C A 7 i 




I N T J J 3 :1 <J 


IF INCS**. t„'.l I uC TC 501 




IN TCo -SO 


fOC REaJ ( L, l iz\ r p A 




IN T J JJ 


WK IT £ (o , 1 jo ) 




1 N TO 09 10 


PE AJ (5, 1.15) PP 3 




I N i J J J ^ J 


WR ITE< c» lot ) P°i 




INTOJSJJ 


FcAl (5 • 1 io i PFC 




1 N T C C S 4 0 


kR 1 TE( c , l J .) PFC 




1 NT G o 0 b 0 


PEYC (5, leS ) ,0 C 




IN7 G 09 c;0 


URITL(o,l.io) PP„ 




IN TCC97G 


50 1 RE aO ( 5, 1 3; ) PSM 




in r j j j j j 


RR1TE(o*Iao) P^lJ 




IN T0G9S J 


IF ( ISCaA.cv. 1 ) ue T C 131 




in re icco 


18 C IF ( IiC V) .0 1 . 1000) C J TO 30 5 




I NT J u I J 


ir=i 




I N T 0 i C 2 0 


ICDPS-1 




in r j io 20 


ISCA.N=ISCX -it-1 




I NTOiOTO 


FS i i\-P$ i M5 .0 




: INTCICeO 


181 PS^ih 




INTO 1070 


F SC = F S I x 




IN TC 1 C HO 


psr=ps 




I NT 01 OSJ 


C C A LCUL A T t ^ RCESS’JRE COEFFICIENTS 




IN TO i 10o 


30C PJYH=PT-PS 




INTO I I i J 


CPi= ( PP 4 - P j ) / P 0 Y !>! 




IN TO i 1 £ 0 


CF£ = 1 FF8-F5 l/P^YL 




I N T C I 1 J 0 


CPC= ( PPC-Po) /p Or t\ 




in roil 4 j 


C F C= ( F J C - 3 j ) / J CY N 




IN 7ui IbO 


C CGFFECTS Pi ;5 GMPTaCN Tr Efi SORE CF ARE IN CALIdRATlJ.4 RAN^E 


IN TCI ICO 


IFIFOC lPi.c j.2) w J TO 301 




INTO 1 i 7 J 


Z AX' AX = - 1 OOu J . 




i N TOLL 30 


CAlL X'tXA-'YlZAjZAX'AX , N X , NY ) 




IN TO I 1 y J 


IF ( CP A . C f .<_ •••« •* A a ) PaC = 1. 0 l-^S 




I NT J i 0 00 


IF (lFC. J C >C=1 .01 APS 




I N T C 1 c 1 0 


lFICPO.o 1. CAMrtX) Pa^-l.Jl-PS 




INTCU20 


Z E - '<Aa=- 1 JyjJ . 




1 N Tu I c. JO 


CALL F A X J - Y t Y c t A E v AX , IX , NY ) 




IN TC 1^9C 


IFIC^.JI.ZIUX) PS 3 = 1. 01* PS 




1 NTO 1<L 30 


2 A V 1 F = 10JJ0 . 




IN lul^oJ 


CALL FlN4rY(ZM,Z>-FiN,\X, NY ) 




INTO LI ; 70 


IF(CP*.t-T.Z4-ilN) ->j[ = L.YV*»S 




INTOIx_dJ 


IFICFC.Lr.AJ-'iN) F 7C-0. PS 


• 


I N T C I 2 c . C 


IF (CPU.eT.ZA'l t) F 3C = C. yS>*d$ 




IN TO I £00 


Z8F r = i jjjo . 




INTCiJiJ 


CALL f'lLJM (n- , ZE* tN .NX , i\Y ) 




IN TO 12 20 


IFICPP.LI.Z0 '1 1) p-SC ^C.VO^P 3 




INTOlJ Jy 


LF IDSC.CO • J S 1 j C ' 7~ OJi 




1 N TO i J h ) 


ICOP S=1 CwPSf l 




1NT J IJ^J 


IF ( lCHPS.ST.lO) CO TO 180 




I NT 01 JoO 


P S - PS c 




iN TC 1 2 iO 


P S T = P ;C 




I N V l 1 j 1 0 


CO TO jJO 




INTul jNJ 


C CALCLLAlhS T *.} -lF C FHI ANjLES FOF ’)» t 'll’ 


* 


INTO 1400 
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F I LE : VDP 



FORTRAN PI 



p 



S 



N 



3 C i CALL I ^ f s: S ( 7P\ ,YP*,YHLC,Y PUP, WIN, 

1C FA ,C Rd, ALh , ^,U , a A f Y A ,Z A , 

2Xt) »Y t f i E ♦ A A X t Y ^ A , 7 A X f X t X » Y P A , Z EX , N X , MY , IE PS ) 

IL PS 1=15 PS 
AL i 

ALF2 = AlI- i > ) 

Ph 11=PH I ( n 
PH i 2 = Fn 1 i / i 

C CA lC L i_ A 7E 5 Ia. hlF £ PU ANGLES FDk • 1 * E •III’ 

CALL I JTSCj ( XP A, A 3 :, / r L *\ aPU°, X<1 A , 

1CPA , C P L »f.il , A L h ,a U ,Y AX ,Z U , 

2XAXfTAA,ZA/,A4f Y *-» » <i A f / A , Y^-. » Z A » !\Y , NX tit PS 1 
It PS2= I E ° r > 

ALF if ( 1 ) 

AL F 4 = AL F ( c) 

FHI3*F,il (i) 

PH I 4=PH M2) 

C StLECTS The. P-MPtP AlF E Phi TuGLE.5 n-JT OF THE FjUK VALLES 

F A Lh^ l l A L F 1 * ~L r2 + )- A M A.< i ( AL F 1 , AlF 2 , AL F J, 

1*LF 4 )-* 4 1 J M AlM. , t-i lF > ,A Lr M ) 7 2 

fPM=((- ) Mi + , hL^hL^hlt)-VMAi(PHl j,FHH) 

l-AMlnll PH II ,PHi2 , PHI j t 3u l t ) ) / ' 

C CALCULATES JF_.s PS K ! v U AT A Uh PIS IT ION IV 
AlFCMF Hr-X PU 

CALL I M FL U Am, r ^ ,Za ,ALFD T ,FPH 1 , NX , NY ,C PU) 

/L F CN= F ALF-A 

CALI lMFLTlXMtYAtZAiALFCO.FPh J , Na , NY, CPC) 

PSN = (PPG <0 ^0-PP L>*C^C )/ l CPJ-C PC ) 

C CCNVERllNLc TlCTS 

EPbPS^iSl l Fi'^-rS) / p S J 

EPSPSC=.cjuCju^uLC.s^{PT«-2.o-.ooo)oo^a+pr 
ih ( IlFSI • l U • 2 ) lP SP’i ^ = .JO) JO 00 JUJ U*l PT **2. J UOJ JO Jc^T 

IF ( i EPS 2 .£ 0 .*: ) £ P j d S G - • l lJ COOOl Cu 9* ( PT** 2. C) U0JG006* D T 
IF (EPSPS .cl • t PS PJu) LC TO 310 
GO TC 3 C 3 
310 1 IT= I IT ♦] 

IF (I 1 7.0 1. I J ) ^ TO 130 
P3=PS+4ELaP 6*( PoN-^S) 

FSC= PS 
PST-FS 
go rn jju 

303 kKlTL lot I7i ) FALf 

171 FjkMTMCS' VtL CITY YAW IS* ,F15. ?_) 
kRl T E (6 ,!/._) hPni 

172 FUFVATI* WLCuITV PITCn IS 1 ,P15.2i 

kR 1 f c. I o , 1 7 j ) P o' 1 

173 FC^^i' 5 T AT I C FStSGU XL IS 1 ,F15.3) 
kR 1 TE l o , 1 7 4) Pf 

1 7 A F C P v A T ( * i’JTAL P-E^OL’E IS* ,F13 .j) 

AP = SC FT (( l kT / PSA j **J .2d3 7i 4-1 .0 ) <-5 .0 ) 

WRIlEU, lb) A 4 

175 FGF^Tl* UCh 'JJS«ER IS 1 ,F15.5) 

CT CALL I ALL An U ) 

CT V 1 RT I M = ( A (3 ) - 1 Ao )/ / 6boj 

CT CPI - l A ( o ) - iAC i. ) / 7C30C 

CT *R ITE( c, A5^ ) VI % li 1 

450 FO^N/TC Vl.TJ.L T I M E IS' ,F15.3) 

CT kR 1 7E l 5 ,45 1 ) CPU 

451 fJRMATl* M7\L CPU IS* ,£15.5) 

CT iA4C = AiM 

CT 1 A 6 C =m ( 6 ) 

CO T ^ lv 2 

3 C 5 «U TCI 3 ,23 7) 

30 7 fOV-*AT(» ^40 L -C 1000 TIDES' J 
oO TC 10 3 
ENi) 

SO t ? OUT HE MTS Co l \R £> 1, \^^,A<LT, \PUP ,A21M , 

1 H T 1 , H T 2 ,-c^l , ^ ^ S c. ,^1 f V 1 , Z I , 

2X2, Y i,<- t ,/J, 7 3 ,i_ ^ t'T'tti.o t f\ 1 , N 2 t I EPi) 

Cl^tNSlM X [ 5 i ( i J ) , L 3 <. I i J ) , ^ lO ( i j ) , A 1 ( 4 j J f r 1( J ) , 



IN T J Is 10 
1MTC1H2J 
IN J C M 20 
1 iN T U M ^ J 
10701^30 

into r,eo 
INTO 1 4 70 
I N T C M E 0 
I*j f C i 4 SO 

l r , *r o i30J 

IN T C 13 10 
INTO L 320 
IN T U 15 30 
I N T 0 1 *J 4 J 
INTO 15 30 
IN T l 1 5 t 0 
1NT015 n 
IMT0I33J 
IN TC 15 SO 
IN VO MOO 
IN7G161J 
I V J TO I c 2 0 
1 NT 0163 J 
INTO 1640 
1NT C i c 50 
INTO lb tO 
1 N T C 1 5 7 C 
I NT 0 1 C 30 
I N T C 1 c S U 
INI Cl 7 CO 
I N T C 1 7 1 J 
IN T C 1 7 2 0 
I NT01 7j» ) 
I NIC 17 AO 
1NTC17SJ 
I NTO 1760 
INTC17 7C 
1 N T 0 1 7 o 0 
I N TO 17 NO 
IN T C IdOU 
I NT 0 Id 10 
IN T Cl t20 
INT 01030 
1NTC1EA0 
1 N T C 1 6 
INTO L Sou 
1NTC1E70 
IN t 0 1 o 6 0 
IN TOIL NO 
IN TC iSCU 
INTO 15 10 
1 MG192J 
INT CISj j 
INTC iSAJ 
INTCISLu 
INTO iSoO 
I N t CIS 70 
INTO IS dO 
1 NT 0 1 S 90 

rnc2i.cc 

INT C2C iu 
I N T 02 o 20 

INT 0 20 30 

I N T 0 2 u 4 o 
INTC2G50 
INT 0 20 60 
IN T020 70 
IN TC2 J JO 
INT 020 S j 
I NTOc 1CJ 
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FILE: 



r jr f r a 4 n 



P 



S 



1Z K 4C #4 3J ,X2( 4 0) f Y.2 i 404 f Z2 <4 0 #40) «X3 (40) «Yi <40 ) « 

2 2-> (sOf i J J I Ai( 4 J i f ¥ 4 ( tC4fZ‘tH40,4\j) 

Pt a m = > k> j • 

RtSl( 2 ) = 40o0. 

FE^2 U 4 = 9 ) 4 J . 

RE 5^(2) = S 0 10. 

IEP$= 1 
IS 1 = 1 
I S p - I 
CAR = 1 . 

AR2^ = AR i K 
GO 10 152 

150 AP2 r ^ = 7 k2R *lAR 

AR2R J=4X2R 

152 AR M 1=AR2R-AR D 1 
A.\P2 = AR2*-pRP*l 

C CHECk FCh Lirtcr* CALI KANjI 

IF ( a .L T • A'v L ) ) L ) TO 15 0 
IF < A<«F2 • IT.AkLC) l! T 0 130 
K= I 

C CALCULATES IN 1 T I *-L PE-JETRATICN POINTS 

CAtL rcNFTSUl ,.\^ ,<1 ,Y1 ,Zl ,H T l ,ARn , AR*ES J 
RAJL = A RRfcS ( 1 ) 

I'f't ) {2 j 

CALL PEhPTS(..l f i*2 »a 2 »Y2»Z2fHT2 ,ARM2 , ARRET. ) 

P3JL = A^tS ( I ) 

RB j R = AkR t> i2 ) 

C CALCULATE^ SlCCeblVt °ENETRAT I CN PCI NTS 
CO 140 1=1,10000 

AR2R = AA2r< + JA* 

ARM 1 = ZR 2R- AR P 1 
ARN2= AR2R-AKP2 

C CHECK F C p LPPEk u^LaLRmTICK RAKuE 
IF ( A3M.3Y.ARUP i GO VO I 95 
IF (AhN2.o"*.A,-UFJ U TO 1 Sr 3 

CALL PE \ P T C ( \ 1 » N ^ i a 1 f Y 1 t Z 1 # n T1 »AKM1 , A R R d S ) 

RAlL=APLtS l I J 
FAih=ARRt$(2 i 

CALL Pc'IP i j< 41 ^ ,X2 ,Y2 ,Z 2 ,HT2 ,Akh2 # ARFE S) 

PBIL =A^R£S ( 1 4 
R81R-ARR ES( 2) 

I F ( R A J l . L J . iOOo. ) I S L = 10 
IF ( R t J L . tC. 10 J ) • 4 i r ' 1 = i o 
IF(k8js.EC.is>CJ. 4 I 3 R = 1 u 
IF ( R A J R . tO . ^ = iu 

IF < RAI L . EC . 1 ) j 0 . ) i ^ L = 1 ) 

I F ( R F IL .E ) . I J J 0 • ) I SL= 10 
IF(FEIR*t4.i)JJ«) i S R = i J 
IFIRAIk.EC.I 400. 4 I 5^ = 10 

C CHECK FOR JlRirCT HIT 04 LcFT INTER SEC T I ON ,E VA L UA TE IF RECUIRED 
IFIISl.Kt.iJ) „C T J 164 
I S L = i 
GC TC 163 

164 IF IRAIL.Kl.R6iL) lC 160 
RE SI (K)=R 3 IL 
PES2 (K 4 = A R 2 p 
IF (K. £0.2) ^ J TJ 155 

K= 2 

GO TC 13 J 

c check fjr i jte> vfji ate left i ntersect ion, evaluate if required 

180 IF ( laL . 4t. iu) GJ T J it J 
IS 1 = 1 

GO TO If 3 

160 IF ( ( ~AJl-R u JL ) / ( : A I L S I L ) • G T • C • 0 ) GO TO 16-> 

AA=(lr2rJ— / 4^AJL“RAiLi 

BA=AR 2 <- 1A*-. < Ul 
A B= ( Ah2FJ-APJ-> )/ < * IL ) 

B3=AR2^-^ 3 i L 

RE S 1 ( K ) = ( 0 * V J / ( A^-AR) 

RE 52 (K) = A^*RESl(K)i-BA / 



INI C 2 l l 0 

INTQ2 120 
IN TL2 UO 
1NT 02140 
IN Tu> 1 3 J 
IN T C2 1 1 0 
1 : 4 T 0 2 1 7J 
IN TC21 FO 
IN 102 1 fO 
I N T 022 00 
INTC22 10 
I N T j 2 2 2 4 
I N T 02 2 jO 
I N T J 2 2 4 J 
I NT 02 2 30 
IN Tl 2260 
INTCc2 70 
IN TO 2^ 30 
INTC22SC 
1 NT 02 3 00 
i mt: 23 n 
I N T C 2 2 2 0 
INToZi JJ 
I N T C 2 3 4 C 
IN T 02 3 :>0 
I NT 02 3 90 
IN T l 2 3 ?0 
I N T 0 2 3 00 
I N T 02 J 90 
INT 0 2 *i 0 0 
I NT024 10 
4 N T C 2 4 2 0 
INT 02430 
I NT 02 4 40 
INTC245C 
l.N T C 2 H 6 0 
IN Tu24 / J 
IN TO 2 4 60 
1 N T 02 *t90 
I N T 0 2 5 0 0 
INTO ^5 IJ 
I NT 02 52 0 
IN TC2 5 30 
INTC25hU 
I N T 0 2 3 3 0 
IN T 02 56o 
I NT 02 3 70 
IN TO 25 30 
INT C 2 3 5 0 
I NT 026 00 
INTC261C 
INT 02620 
ItN T026 j>0 
IN T02L40 
1NTU2650 
INT02660 
INT 0 26 /O 
INT 02600 
1 N T 0 2 o 9 0 
INT027C0 
INTJ2710 
IN TC2 72 J 
IN V 0 * 7 30 
INTJ^ 74J 
IN TC2 7 50 
IN TO 2 76 J 
INT027 7J 
INTO ^ (80 
I NT 02 7 90 
IN TC28C0 
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> 



FILE: VDR Fu RfRAN °1 



N P 



S 



I F < K . E 0 . 2 ) T J 

K= 2 

C CHbC* E LR ui RtCT HIT 

183 IF U jP • 7b. I J) ut) 
13 «= 1 
GJ TC Ki 

U 1 Ir< RA K . IF K UR ) 
R E r > 1 U )=R8 1R 
«E S2 <M = A -<2 < 

JF( K .t J .z ) u-) T J 
K = 2 



135 

C) RIGHT INT fcRS ECT I)N, evaljate 
TJ loi 



GJ TJ 102 
iot> 



C CHECK FJR [\\6V4cDl\TE R IGHT l N TE R SEC TI UN ,E VA L UA TE 
182 IF ( I iR. At. 1 J) oJ TJ 162 
I S R = 1 
GC TC ifej 

162 IF ((PpJP-HcJO/( PM^-^EIP ) .GT.Q.J ) G C T"i 163 
AA= { AF2RJ-AR2P ) / i K ) 

fc A = 



A3 = IAF2' ? J-AR2k)/ ( KjR-Rol -l ) 
PB= AR2<-A-R«U-« h 
Fcsi i k i =(Fi-eA)/ 

RE S2< h ) ^Ah^uoI ( iO f 6 3 
IFU.EQ.K C- ^ TJ 15 5 
H = 2 

C RE-INI VI ATE 
163 A82fJ=A^R 

R A jL =F. A iL 

PAJR = F<a1R 



IF 



IF 



REO'UREO 



RtCUI RED 



pe.;l = r d IL 

14C R8JK-P31R 

fcR ITEIo ,2iJ ) 

210 F 0 r ^ AT ( 9 lJCP 14 j THROUGH, NO PCINTS FOUN C • ) 

GJ TC 156 

195 IF { RtS 1 ( J I .c 0 . 9UCC. U G'J TJ 191 
IF IntSl (2) .cC. >000.0 ) G c r C 191 
I H R t f, 2 ( N.cJ.SClC.C) 1J TO 19 1 
lF(RtS2 U) . E0 .9GJJ.0 ) G 1 TJ i9 1 
GO TL 155 

C EVALUATION JF MTERSbCTIJNS *HtH CLRVtS ARE ALMOST TANGEM 
191 YGOO.O 
I E P S- 2 

CALL P£'!P7g<I2,N1,Aj,Y;>,Z 3 , HT 1 f YGH , A PRE $ ) 

ri = aes ls 12 )- a^cKI ) )/z 

CALL Pt JP|C( J2tNl ,Y^.Z4 f riT2»YGij f ARPES) 

R2 = t 6F ( AR US U )-ARR ES ( l) ) / c 

RE S2 (1) = (( RK* 2) -( R2 -<2 ) + I ARP2**2) ) / (2.0 * ARP 2 > 

RE32(2 ) = R r G 2 ( 1) 

Pb SI ili =3Gr<T ( AeS ( ( Ri **2 ) -( -ES2 (1 )**2 ) ) ) 

RE S 1 ( 2) *-*l SKI) 

155 RETURN 
E N J 

SU3R 1UTI 4E PENP T S ( IX, NY,X , Y , Z ,HT f Y G , X R ) 

01 v L ^ S i J N i KJ J , Zx( 4 J ) ,Z ( w, VO ), XI A j ) ,XR( lJ) 

18 = 1 

C SEARCH fop J JF lJ.xEn i LINE 
0 C 101 1=1, \Y 
IF ( Y ( I ) . Sr . Yu) >c rc 102 

101 CCNTIMJc 

102 0=1-1 

C COaPjE A \ F ZERO o\S5 

C NEXT A LTicK ARE EXECUTED IT FIRST SEARCH ONLY 
1C = ( to- Y ( Ji ) / ( Y ( sjtl j -Y ( J) ) 

Za( 1 ) = I Z ( 1 , J* 1 )-Z ( 1, J) ) KC + K 1 , J) 

IF (2X U ) . fct*.hT j CC TC 1J7 
13=2 

112 00 1U 1=13, 4 X 

ZX ii ) = ( L ( I , J ♦ L ) - < (A , J ) J*Y C *Z ( I , J ) 

/^(Z(I f Ju)-Z(i-i,j))M Y,-YI J) ) /{ Y ( J+l ) -Y ( J) j *Z( 1-1 , J) * 
IF IZX ( I) .£0 .n D lJ T J Iuj 



INTC^E u 
I N 702820 
IN TC2 fc 20 
INT02FAJ 
IN 702 3 50 
HT C2ct0 
INT-J2 3 ?0 
lNT02tcC 
INTOKi U 
I NT 0^ 9 00 
INTO 29 10 
HTC2920 
1NTC29 JJ 
IMTC25R0 
INT J2 950 
IN T 02 9 60 
IN T C 2 9 7 0 
INTU29 3J 
I N T 029 90 
IN TOJOOU 
I NT 03 u 10 
IN TCSC20 
INT C3u30 
I N T 030 40 
IN TC3L5C 
INT J J JcO 
IUT03070 
IMTC3CEJ 
INT03J9O 
INTOjiOO 
INT03 1 1J 
I NT 03120 
IN TC3 130 
INTC314u 
I N TC J 1 50 
IN TC3 160 
INT031 70 
INTJ3 180 
INT 03 190 
I NT 03 EuJ 
INTC3210 
INTO 22 20 
INTC32 j 0 
IN T C 32 40 
I M TO./ 2 50 
INT03260 
INTO j 270 
I NT 032 80 
IN TC 22 SO 
INT033C0 
IN'TOiJ 10 
INTCJ320 
IN T 033 30 
IN T 03 340 
INT 03 3 50 
in f 03360 
IN TC33 70 
INTO 2 2 dJ 
I N T 03390 
IN TC24L0 
IMTOJ^ iJ 
intoj42j : 
INT C 3 ^ 3 J 
I NT 0 

INI C 'T50 
INTCJ^tcU 
INT 0 J^7J 
IN TC34 3 C 
I N T J 3 a 9 0 
I N 7 03 5 JO » 
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oro»— 



File: VJP 



FORTRAN PI 



N 



P 



S 



IF C ( ZX( I-n-iHT) /IZXUI-riT) .LTbOvO) oC TT 104 
10 J CONTINUE 

vj C T C 1 1 o 

C CALCULATE JF »ENG |'F ATI CNS IN fc ITIIER TRIANGLE 

104 x^xii-iim^i n-xii-i)i*(ifG-ti jii/ivijai-nj) i 

IflZR.tO.hT ) U C 1C 105 

If ( ( ZXI I )-H I) /( Z.i-HT ). uT.0.0) GC 7 C IC5 
LU = A t c (f T - t ZX (l-i)-HT ) 

>£=)(! - I ) + ( x ( i — 1 > ) <“ AbS I ZX ( I -I ) -HT )/UD 
GC TC 110 

105 XS=X V + iX i i ) -aN)* AdS ( ZV -hT )/ ( APS IZ'-ET M A85I ZXI II- HT )) 
G j in lit 

C LOGICAL FtNETA ATIUNS ACC ONUL A f ION 



INTC3510 
INTO 3520 
I M T C .1 5 3 C 
INI J 3540 
I M TO 3 5 50 
INTOBdcO 
INTJ357J 
IN TO 3 5 60 
IN TO 3 5 50 
I N t 03o00 
IN (0 3 6 10 
INTO 3b2 > 
IN TO Of 30 



1 C 7 


>K(1)=A|1J 






IN TC3t4C 




1= 1 






19103050 




GC TO ilJ 






IN TC3t.6J 


lot 


4 1 R ) =- x C i) 






I N T 0 3 1 7 b 




GO 10 111 






I N TO 3t 30 


110 


X R l I R J- X S 






INT0J690 


111 


IFI1R.B7.2) ^0 


ro 


114 


INT 03 7 00 


113 


1R = 2 






IN TOb 7 10 




ie=i+i 






INT C 3 720 




GO (D 11 < 






INT03730 


116 


IMIR .tC.l ) GC 


TC 


ill 


INT C3740 




XR C2) = 10JJ.0 






INT 0 3 7 50 




GO TC ih 






I N T 05 l bO 


117 


XR (1J = 1 JJJ .0 






I I TC3 7 7C 




X<i ^)- 1 l L C . G 






IN 605730 


1 1 A 


RETURN 






I N TO 37 90 




END 






I NT 0 3300 



SU"*ROL T I Jt I 9 fPLI ( X ,Y , Z ,Xb ,YG ,t\X,NY, ZRES) 

C LINEAR iNTSRFCLA T ION CN CAL I BR AT I IN 5 U 0 FACE TO EVALUATE 
C CP VALLE a T Ct-uKUl N aTcS XG.Yb.THt PESULT RETURNtO iS ZRES 
01 PE NS ION A<4d),Y(4C)»Z(4b,4C) 

C SEARCH ECR 1 CE LcET X LINE 
00 1 1 1 = I » -4 X 

IFIxIIO.CE.aGJ GO TO 2 
CCNTIME 
1= IC- l 

StAPCf TCP J CE LCWE? Y LINE 
00 3 JC = I * j Y 
If ( Y( JC ) .Gt .YG ) GO TO 4 

3 CONTINUE 

4 J= JC - 1 

AL= (Y (J ♦ ! )-Y ( J I )/( XI i + 1 )-X( I)) 

B L= Y IJ) -Jl-X ( 1 ) 

YCR=(AL * > G ) ML 
IF (YG.LC.YCR ) GO T) 7 
IF ! Yo . •> T • YC -\ ) to Tt 5 
C RESULTING Z IS II L JRRER TRIANGLE 
XJ=X ( 1 f 1) 

Y 3= Y ( J } 

23= Z 1 1 ♦ 1 1 J 1 
GO TC 6 

C RASUL f INC Z IS IN LF> P ER TRIANGLE 

5 X3 = X ( 1 ) 

Y3=Y( JU ) 

Z3=z ( i, j + n 

C GENERA l CALCULAT 1 C M : tCPC FCR eCTE TRIANGLES 

6 xi=xm 

Y i=Y ( J I 
Z1 =Z (1 , J) 

X2=X( I 1) 

Y2 = Y < J ♦ 1 I 
Z2 = Z< ML .J»ll 

AP= ( ( Z l— Z J ) * ( Yl-Y2)-(Zl-Z2)*< Y l-Y 3 » 1 / ( ( XI -a3 ) * ( Y l- r 2 J - I X l -X2 J 
1*<Y1-Y3 ) ) 

6P = ! ( ZI-Z2I- \i>* ( XI-X2) ) / (YI-Y2 ) 

CP = Z1-AP*a i- J .P*Y L 

2RES = AP*X u *■ EE* Y t *Cr> / 



IN TO Jb lb 
INTCJS20 
IN f 0 3 8 3 0 
I N T C 3d 40 
1 INTC365C 
IN T 03 360 
l N TOO d 70 
IN TC 3860 
INTO 3050 
I NTC35CC 
I NT 0 39 10 
I NTCjGZO 
I NT CoS 30 
I9TC3S40 
IN TO 39 50 
IN TC 39 bO 
1 NT 03 5 70 
IN TO 3930 
INT C3S SO 
I NT 04 GOO 
1MTC4CIC 
IN T G4020 
1 N T 0 4 0 3 0 
INTC4C40 
I N T 0 4 0 3 0 
I N T 04060 
I N T 0 40 70 
I NT 04 0 80 
IN TC4C50 
IN1C4 ICO 
INT041 10 
I N T C 4 I 2 C 
IN f 04 l 30 
I N 7041 40 
IN V C 4 I 50 
I N i' 04 1 b J 
INT041 7 3 
IN T 04 IcO 
I NT 04 1 bO 
1NTC42CJ 
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Flit: Vo« FORTkA J °i N P S 

GJ 7 C LOC 

C RESULT l 4 Vi L IS M JIVIOING Lite 

7 xres = u i in , jh j -i 1 1 , j 1 1 mxu-x i m/ ( xi i +i >-xi m *z ( i , ji 

IOC c E f URN 
ENC 

S L cR J L I I N L R yXARY I Z , 2 W AX , NX , NY ) 

C EVAL JA T ION 7F M*Hu< CP Jt- CALIBRATION SURFACE RaTRIX 

twf' ii-tOjAiJ) 

JO J I = 1 , NX 
00 i J^i, IY 

i If IZiJ , JI.oE. Z^AA J Z w AX=Z(ItJ) 

RE TURN 
ENC 

SLtROLTIiNE "INmSY (Z , Z R I IS , NX , NY ) 

C EVALUATION OF N I n| IN J CP Zt C a Ll 6k ATI ON SURFACfc 1ATR1X 
CIfE>Sl.N tltorul 
OJ I I - 1 t it x 
CO 1 J= 1, l> 

1 IF IZ I 1 , JI . Lc. z^l.\ ) Z»i J=ZI I, J I 

RE turn 
ENC 

C END Cr VCR 



IN T 0 A 2 10 
INTO A2 20 
IN UA2 iC 
IN T0A2 AO 
1 N 7 0*t 2 l> 0 
IN T C A 2 1 0 
INTO A 273 
I.NTCA2EC 
I N 10 a 2 '> .3 
INTCtJOO 
IN 7CA3 10 
INTO ,3 20 
INTO yJ iO 
INTCy3A3 
I NT 34 j JO 
INTCA JCO 

In i c so ?o 

I N i j A 2 J J 
IN T C A3 SC 
IN T 0 *A0 3 
IN 7 0 A A 13 
IN T C A A 20 
INTO ♦ A 30 
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APPENDIX 2 



LISTING OF SWVDR 



U/Ct//5 it. 45. IS 



Flit: S-mVCR 



F CRT RAN Pi. 



C 5 


TAR I 3F SaVJR 




INTCuC 10 




Ci^ENS 1 H ->LF ( 10), XA ( 4 ( ) , YA ( 40 , 




INTO JJ^O 




1 L A 140 .40) , <d (4 0 ) , 1 J (4 J 4, L HtO, -MJ ) ,X AXi 4J ). TAXI id) 




INTCCC2C 




2 Z A > ( iC) » x a ( , 1 3 x { h u ) ,Z.3X(4 0,40) 




IN T J ) J/,0 


119 


FCR* 4 AT 12 15 ) 




I N TOO 030 


131 


EursVA TC f 7 • 1 J 




IN VC 


131 


f-TV-IAUf 5. 1) 




I • 1 T J 0 J 7 0 


120 


EC h v ( E 6 .3 ) 




INTCCChO 


132 


ElH^A 1 1 F I.), z) 




INTO JJNJ 


135 


FCP‘*1AT( FitJ.D) 




1 NTJjIOO 


136 


PO.* AT ( E l 7 . 5 ) 




IN 700 1 10 


C HEAD P v J-<t C M. 1 =U ATI .J 4 DATA 




IMT JO l 20 




REZC(2»il-y) > 1 X , J r 




INTOJl 40 




ON I 2 i 1 =1 • JX 




INYO J 1 TO 




R E 4 Q ( 2 , IJj) XA Cl) 




I iJT J Jl^J 


121 


C',M IiVJE 




INI COl tO 




DJ12^l=l,fw 




I N i o C 1 7 0 




RfcAOUt lij ) YA( I ) 




in rooi <n 


122 


CC M i <\ j t 




IN TC Cl 50 




OJ 122 1=1 f NX 




IN J2l0 




CO 123 J = i, U 




IN 1 0 0 2 1 J 




REZO (2,120) ZA( 1 , j J 




IN T C C 2 20 


122 


cow i inie 




1 1 J T J Z 4 J 




CO 221 1=1, MX 




INT0J2 AO 




RE AJ ( 2 , 120) Xri ( I ) 




IN TO J 2 50 


22 1 


CCNf IMJE 




I NT0J26O 




CG 222 1=1 , NY 




IN TL02 iJ 




RE AD ( 2, 120 Y6( 1 ) 




IN T 002 60 


222 


CJNT IflJE 




I N T 0 u 2 5 J 




DC 223 1 =1 , NX 




INTCC3CC 




03 2 ci J=i, Jr 




IN T 003 10 




READ U, 120 i L E( I , J ) 




IN T0032U 


223 


cun time 




IN TOO J jO 


C TRAN6EGRV CALInR^TlJu MATRICES 




INT0u34j 




cn zz 4 i=i,ny 




• 1NTCC35C 


224 


>a x( i ) - y~*i i ) 




INTOJEcJ 




00 t25 1 = 1, Ia 




I NT 003 70 


225 


VAX ( I )=X 4( 1 1 




I N T 00 2 60 




OJ 22b 1=1, \Y 




1N7CC3SJ 




on tit j = j,'u 




INTOUhJJ 


226 


ZAXli , J) =Zp1J,1 1 




INTCCA1C 




JJ 227 I = 1 ,,\ 1 




INT0GA2 J 


227 


XdX ( 1 )= Tt (1) 




IN TO Oh 3 J 




00 22b i =1 ,NX 




IN T CCA AO 


226 


Yd > ( 1 ) = X P ( 1 ) 




INT0G45J 




CO 225 1=1 , TY 




I N T 00 A 60 




00 229 J = 1 , MX 




INTOOA 70 


225 


Zixl 1, J J=Z°( J.I) 




I NT Ou A 60 


C FRCBE SfcTTiMGS 




IN TCuA 90 




YP A = C. 0 




INT CO 5 CO 




YP 0= 55 . 




INTOoblO 




X P A = 0 . 0 




I.NTC0520 




XPC= 20. 




INTO J 3 20 




XPC=-2J . 




I \ TO Jb AO 


C PKCdFAF CC.NSTfl.MTS 




H 7CC5 50 




L0=-JO 




I N T 0 u 6 o u 




Y RU P = 3 J • 




I N TOO 5 f 'J 




>k IN=-3 0 • 




I N T 0 v> b 8 J 




XRL0--3 J . 


* 


I NT 00550 




X R L P = 3 0 . 




IN T C060C 




>*IN=-3C* 




INTCU610 




1C0PS= 1 




1 N T Ojo 20 




NCC C FS =2 




IN T C C t 3 C 




TJQSI M=2 




IN TO >6 AO 




1 1T= 1 




IN T0oo50 




1 SCAN=1 




INTCucEO 




REl>P6=C.3 




INTOJo/O 


C EX P E R 1 M E * J T 3 1 1 JL AT 17 M 




IN 7006 30 


153 


C'IN ll^lc 




I N V OOoOO 


401 


f JP ^ AT ( E i 5 . h ) 




1 NT 00 7 00 
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FILE : SWVDK F CR1 r«N PI 

REA0IS,4Cl) FT 
WR IT E ( c » *,0 l ) PT 
IF (KSi C.uJ ) oC TO 500 
PE AO ( 5 1 4 l L) PS 
WRITE 16 » -tJii P$ 

Rt<OI5,«CU \LI-A 
wRlTc!c*4(,l) al F 0 
PEAC (5 ,401 j fHI 
wRl Tt ( i ,401) Prli i 
CALL I'JTPL I I a A, Y A, Z A, ALFA ,PHII ,NX,NY,CPA) 
ALFC = ALF A-a PC 

CALL 1 1\ I P L T ( X 4 , t ,, ,2 a ,ALFC,PHI1 , NX, NY, CPC) 
AL FL= AL M-XPC 

CALL I MFL7 (a A ,Y 7 ,ZA ,aLFC, ^F I I,NX, NY, CPC) 
P F 1 0 =P H 1 1- YPT 

Call I NT pl t UH,Y F,ZB, ALFA,PHia, NX, NY,C°3) 
PO YN=PT-F3 



N P S 



P P A= C P 4 * P J YN + P S 
FF6 = CFs=*f.,YN+PS 
PPC = CPC* P JY 'ifPS 
F F C=C PO H P CY I 

C READ MEASUhcFfcNTS CATS 

IF (MO SI H.E ). 1) > I TO 501 

500 PEAL(5,i35j p=A 
WR1TE(4,HM PPA 
READ! 5, Ljj) PP) 

WRITE (S ,Ijo ) PP3 
READ ( 5 , I i 5 ) ?PC 
WRITEI j, Ijo) PPG 
RE AO (5 , 1 35 ) P PC 
WRITE! c, I Jo) PPu 

501 FtAC(5,iJs) PS IN 
WRITE! 6, 136) PSn\ 

IF ( ISC AN .GO . I ) jO TO lal 
180 IFliSCAF.^T .lOOO t TO 305 
1 1 1= I 
ICCPS=i 
I SCANMSCAN + l 
PSI N= Pil NF 5. 0 
WRIT E ( 6 * lsfc) P3I.J 
1 e 1 CCMIMt 
PS=PSI 'i 
FSC= PS 1 N 
PST=PS 

C CALCULUtS PRESSURE COEFFICIENTS 
300 CONTINUE 

PO YN=Pl-PS 
CFA= (P?A-PS)/PCYf 
CFE= 1FPE-FS) /PCYlv 
CPC=(P?C-PS) /'>OYN 
CPL= ( PP C-°S ) / p CY N 

C CORRECT' P 5 ASSLPPTIlN TC ENSURE CF ARE 
IF ! NOCOP S.tO .2 ) SO TO JO 1 
2 AF aX=-1 JOoO . 

CALL MA XiriYI ZA»ZAF4X , NX , NY ) 

IF iCPA.CV.ZAMA) PSC=i.J i*PS 
IF (CFC. u T. it -AX ) F S C = 1 .0 l « FS 
IF ( OP 0 .ST .2 \ HA X ) p .SC - i.O i-*FS 
26 Pf X=- 1 J )00 . 

CALL FAXArY 1 26 ,ZEN AX ,NX ,NY ) 

IF (CPE.GT.23 Ux) =>:>C=1.(H*PS 

2AFI N=i)OoO. 

CALL MI N A r> Y ( 2 4 » 2 A M i, * N A* NY ) 

IF (CPA.L T .2 A 1 IN ) °JC-0.99*PS 
IF (CFC.LT.ZJ’MN) PS C =0.5 9* PS 
IFICPD.Ll.ZAHI I) PSC-C.GN-'PS 
2 6* IN= i ) 000 . 

call MINARY! 26 ,ZoMN ,N< , NY ) 

IF ( CPF «L T .Z3N IN) p SC = 0.99*PS 
IF (PSC. EL. PST ) lC TC 3 Oi 



INTCC7 10 
INTOJ 720 
IN TCl 73C 
IN T JO 740 
IN TJJ/50 
IN TC 0760 
1NT0S7 70 
IMTCC70C 
IN !' JO 7 ^O 
I NTCGJO ) 
IN TCGS L ) 
l NTC082J 
I J « J i; T iO 
I I TCCfc-YC 
INF JU 650 
l.NISSobO 
IN I Cue 7G 
I N TO 03 SO 
INTCC8S0 
IN TO 09 00 
I NT 009 10 
IN TC0520 
INTO0930 
1 N T 000 40 
I NT oj 9 50 
I NT 0096 0 
INTCL570 
INTCOGeO 
INT0J990 
IN TC 1 C CO 
IN TO 10 10 
IN TO! 020 
IN T C 1 0 30 
IN TO 10 40 
1NTC1C5C 
INTClOoo 
I NTO10 70 
IN TCIC80 
INTC10 00 
I M T 0 1 1 0 J 
IN TC 1 1 10 
INTO 1 1 20 
IN TO 1 UO 
INTO 11 40 
INTO 1150 
IN TC1I6C 
INT01 1 70 
I N TO 1 1 30 
IN TCI 1 SO 
INT0120) 
I N TO l 2 10 

IN C 7L i 6RAT ION RANGE INTC1220 

I NT 0 1 2 30 
IN TC 1240 
IN TO 12 50 
INT012OJ 
INTC1270 
* INTO 12 BO 

IMTC12S0 
INTO 1 SCO 
INTO IS 10 
IN T Cl 32C 
IN TO 13 30 
I N T 0is40 
INTC1250 
INTOisOO 
1NT01370 
INTO 1 2 SO 
1 NT 01 J 50 
' INTC14C0 
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FILL: 5WVLR FORTRAN pi 



N P 



S 



301 
C ~ 



IC0P$ = IC£PS + 1 
IF ( ICHPo .GT . 10 i TG 180 
FS = P S C 
PST^PSC 
OG TC J JO 
II CGMiMt 

CALCULATE iw) AlF fi PH 1 ViliLcS FC 

CALL 1 K * C3 { ( P A , Y D i f YP L i f y RUP r YK1 n , 

1CPA fCP‘j ,AlP v PH 1 

2XR,YR,Z J . ,Aa f VAX,ZU,A!jXtYHXfZeX,\X f 
ALF1 =ALF ( i J 

A 1 C 3 =A < C I J I 



o'! » 1 * 

I N , 



tAntri'tt » i a ■-* ■ 

ALF1 =ALF l 1 1 
AlF2=AlF (2) 
FhU'PHI 1 ) 
FHl2 = Frll 12 ) 



f. ’ll’ 

NY) 



Kim-iTii ii 

FH 12 = Fell 12 ) 

CALCULA1FS I « J ALF E P, 1 1 ANoLtS 
CALL i NT CCS (XP A ,XP Z , XR L0, API. 
1C PA ,CPC ,PH1 , -cF f , y t X , t AX , 



ALF3=AL+ U ) 

ALF 4=p lF (2) 

FHU=Pm (1) 

PHI 4=PM (2 ) 

C SELECT S THE P<)PfR A or & PHI ANGLES DJT CF THE FOR VALLES 
FALF- ( < U H ♦ -L F t ■* AL F 2 ♦ AL F-, I-AMAX1 ( ALF 1, ALF 2, AlF3, 
1AlF4)-a v I., 1 ULFl ,AlF2,-,LF3 ,ALFh ))/2 
F? FI=< ( PF l 1 +PrlI 2+P 1 1 j + ,; Ml 4) - A 'AXll PH 1 i , PHI 2 ,?‘U 3 ,PhI4 ) 
l-Af'lMlPH 1 ,PFu,Filj,?hl- t ))/2 
C CALCULATE I c * Pi FkH-I OaTA CF POSITION D ANO CORRtCTS PT 
ALF CN= r AL F-X J 0 



CALL 1 J I PoT( Xu ,Yu ,La ,ALFDk,F FH I, NX, NY, CPC) 

AL FCN=FAlF- XPC 

CALL I I T01.T (X A, / A , 7 A , ALFC ‘I, F PH It TX , N Y r CP C) 

PSH = ( P PC *C P J- PP J* CPC) / UFj-CPC) 

907 WR 1TE ( a, 3 I 1)1 IT, FALF, FP,U , PiN 

911 FUPFATl’ 1 T c. i< jT 1 UN VALUES’, Id, 10X,2Fo. 2,F10.2) 
C CfJNVEr.GENCii IE-.TS 



EPSFS-pPS ( ( p SJ-PG I/PS ) 

EP SP S6=.OOOOOOOOCON9 < ( ?T**2.0) - .0 100 0093 *P T 
IF ( ilPji.t 0.2) rPS3S3=.cC^CCC0CLG9*( F7**2. 0) -. C0J>0 003*PT 
IF ( 1 EPS 2 .EC. 2 I t PS 3 G G= .OUJGuJOJUJO* ( PT** 2.0) - . Ov,OOCC£*pT 
IFIEP'PS. jT.CPjPSJ) 3C TC 310 
GO TO jJJ 
310 CONTINUE 
111=11 1+1 

1FI1IT.L-T.ijj jN T 1 180 



FS=F S+r.c.LXPS* IPSH-PS ) 

PSC=PS 
FS T= PS 
GO TC 3 OC 

303 WP ITEIo, 1/1) FALF 

171 FCRFATCO’,’ VlLLCjTY YAW IS’ ,Fl3.2) 

WRITE ( u , 1 /2) F Pn j 

172 FCrMATl’ VlLJCiTY J ITCh IS’ ,F15.2) 

WRlTt(6,l/J) PSN 

173 FORMAT!* i | A i 1 C prcSSO^E IS’ , F 1 5. 3 ) 

WRiTE<o,l7,j PT 

114 FJRMAK* TjT^L PstSSLRE IS 1 ,Flo.3J 

AM= S C FT ( ( l PT/ P Gi- 1 1 **0.285714-1.0*5.0) 

WRITE lo ,1 /p J A* 

175 FORMAT! • MACrt HUK9ER IS’ ,F13.5) 

CO T H 193 

3 C 5 WRITE(o,307) 

30/ -FORMAT! ’ oCwNNEl) UOG TIMES’) 

Lll TC 19 3 
ENO 

SUGRCUT 1 mF I ITSUj! .RP 1, A.RP 2,ARL0, ARUP,A= IN , 

1 H T l , HT 2 , I*lS l , nCjtiXi ,) 1, 21 , 

2X2 , Y i , 2 2 , a j , Y 3 , 2 2 , X 4 , Y 4 , 2 4 , N l , N2 ) 

ClRtNSiO'J -< L G 1 ( 1 j ) , p ES2( 10 ), A PRES! 10) , Xl ( 40) , Y 1 ( 1C) , 

121 140,40) , X2 140) ,Y2 (40 ) ,Z2 !40 ,-tJ ) ,X3 !hO ), Y3( 40) , 



INTO 14 10 
INT0li2O 
IN TC143C 
INTO 14 40 
IN T014 50 
lNlLiiCG 
INTJ147J 
lHTCl-,90 
INI j 1 ‘i ) 0 
INTOIdUO 
IN T C l 5 1 U 
I N T 0 1 3 2 0 
I f) 7 o l 5 30 
I N T C 1 5 4 C 
INTO 1550 
INTO 1500 
INTC157C 
IN TO 13 80 
INTC15S0 
IN V 0 1 c oO 
INTJ16 10 
LNTC 1620 

i NTO l 6 30 

I N TC lo 40 
INTC165J 
I NT 01 tUO 
IN TC U 70 
INTG163J 
INTC1690 
INTO 1 100 
I NT 01 / 1 0 
IN TO 17 20 
INTO 1 720 
INT01/40 
IN TCI 750 
IN TO 1 / jO 
INTC1770 
INF Cl 780 
I NTO 1 / GO 
IMTC1EC0 
INTO 13 10 
I- NTO 1820 
IN TO Id 20 
INTO 1040 
I N T 0 1 o 30 
IN TO IH AO 
INTOia 70 
I NTO 18 30 
INTClcSJ 
l NT 01900 
IN TCIS 10 
IN T 0 1 9 20 
I N T 0 i 4 3 J 
IN T C 15 4 0 
INTJ195J 
1 N T C 1 5 6 C 
IN TO 13 /O 
I NTO 19 80 
INTO 15 50 
IN T 020 JO 
1 N TO 20 10 
I NT 1)20 20 
I NTO 20 JO 
I N TC2040 
IMT02C50 
INT020o0 
IN7L2L7C 
INT0 2CSO 
I N T 02 090 
INTC21C0 
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FILE: SNVCP 



FORTRAN Pi 



N 



P 



$ 



2Z J( 40 fiJ) , X4('i0) f Vf (iC) tZ*(40**6O) 

RE31(ii=*jua. 

Rfc SI I? J =9000. 

Rt32( I) - SC OC . 

FES2 (2 1=9 jO J . 
lSl = l 
IS R= I 
l) A P = I • 

AK2R=A4l N 
GC TC 152 

15C Ak?R-=A*2*+bAs 
AR2RJ= U2? 

152 *K‘U= AR2R-ARP1 

ARM2 = Ar\2 A * A R r2 

C CHECK FOR l'J*c.< CALIBRATION RANGE 
IE ( A c /* t . L T . AhL Gi LC TC 1 5 J 
IF ( AkF2. L i. ARLL) gC IC 150 
K- 1 

C CALCULATES IMTiAL FLKETFATlCN PCIMS 

CALL PLNP r S( U m2 .Xl , Yi , Z WnTl »AR”1 , APRES) 

P A JL= A R R ES ( 1 ) 

RAuR=ARR£S(2) 

CALL P E N P T S ( ’ J 1 , N * fX2tY2tZ2,HT2 ,ARM2 , ARRE SI 
KBwI = /.r RES tl ; 

rb jr=akpes ( 2) 

C CALCULATES SUCCESiVE D EN E TRA T I ON POINTS 
LU 1 AC 1 = 1 ,10 000 
AP *R=AR2P+ J A^ 

ARM1= AK2R-ARP 1 
ARV2 = A*2 K-/^P2 

C CHECK fOP U P*FA CALIBRATION RANGE 
IF ( A FF i . L-T • A "*IJP J n TC 190 
IF( AKF2. ST.A^UF) ^0 TC 195 

CALL PLOPToI a,N2,Xi, r i,Z l,hTi,ARNl , ARRE S) 

RAi L = A P F 1 3 (1 ) 

RAI P = A r< A E S l 2 ) 

CALL PtNPTS ( 11, U, a 2,Y2,Z 2 ,HT2 ,A RM2 f A RRE S) 

FBI L=APP tS (1 ) 

RB 1 R= ARR c S ( 2 ) 

IF CPAJL .E'J.i 000 . ) IS L= 10 
IF IREJL.EG.l ) 0 0 . ) ISL-10 
1 F ( R P J ' < .EQ . 1000. i I S ^ = 10 
IF (RAJA. EC. UJJ . i I C R = 1 J 
1FCP.AIL.EG. 10CC. ) I Sl = iu 
IF( Rt IL .t'O . luo J. i I ^L= 10 
IF (RElR.cC.lO )0. ) i S '- = 10 
IF (RAH. 63. iuGC. J I ^ R = 1 0 

C CHECK FCF DIRECT FIT 71 lEET 1 NT E 3 SEC T I ON i E V4L uA TE IF REQUIRED 
IF(iSL.fsE.lJ) opJ TC 1 64 
I SL= I 

cc tc iao 

1 6 A IF (RAI l.NE.R HIl) GO 1C 180 

RbSl ( K J= A3 * L 
FES2 (K) = AA2 A 
IFIK.EQ.2) G J Tj 155 
K- 2 

GC U 133 

C CHECK FOR l j IEA •itDiAVE LEFT INTER SECT I IN, EVALUATE IF REClIRtU 
180 IF (ISL.NE.iJ) GO T ) lfcJ 
I S L - l 
GO TO 183 

160 IF (( -/JL-PBJL )/( h £ IL-PEIL ). GT .G.Oi GO TO ibJ 
AA=( AR2 PJ-~ \2H / ( A 4 J L-rsA I L ) 

BA^A»c 5 -AT *: .IL 

A8=UF2 JJ-AM*)/ (CEJL-RBIL J 

BB = AK2.<- a 3* k a L 

RES 1 (K i= ( A U ( A v- Vd ) 

RE S2 (K) = A A* F<_$1 IK ) + EA 
IFIK.E0.2) j J TJ 155 
K= 2 



INTO* 1 10 
INI 02U0 
INTC213C 
INT 02 140 
I N T 02 i 50 
INTC2I60 
INI 021 70 
IN TC21 EG 
INTOclSJ 
IN T C22 JO 
INTC22 10 
IN T0L2 20 
IN 7022 jO 
INI u 2 40 
I NT 0 2 2 50 
INTC2260 
INI C22 >0 
I N 7 022 BO 
INTC22SC 
lNTu^JJO 
IN T 02 0 1 J 
INTC2 J2U 
INT0233J 
INTC234C 
IN T 02 3 50 
I N T C2 3 60 
IN TC2 J 70 
IN T 0 2 _> 8 J 
I N T 02 J 50 

rn c 2400 

I NT 0 24 1 0 
I N T 0 2 4 20 
INTC2420 
1NT0244J 
IN T C2450 
INTO 24 tO 
INTC2A70 
IN 7 L 24 80 
IN T U249J 
I NT 02500 
INT025 10 
I NT 02 5 20 
1N7C2520 
INT 02 5 40 
INTC25dO 
INI C2 5 E C 
INI 02 5 70 
I NTC2580 
INTG25SO 
I N T 026 00 
I Nilotic 
IN TO 26 20 
I N 7026 30 
IN T L 2 1 40 
INT 026 50 
INT02660 
IN 7 C2670 
IN T0Z6fcJ 
I N 7 \J<S o 90 
IN T 02 7 oO 
INT027 1 j 
INT C2 720 
INTC2720 
I NT 02 740 
IN7C275Q 
IN TO 2 760 
INTO 2 770 
INTC27EC 
IN T 02 79 j 
I N 702600 
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FILE: 


SaiVJP FORTRAN 3 I N P S 




C CH EC is FO* DIRECT HIT 7l XIGhT I NT fi,*G£t.C T lCN»CVALUAf c IF RtCUlKfcC 


IN T C2o 10 


ia j 


IF ( I SR .Nt . 10 ) GO T 1 16 1 


INTU2623 




I S R = 1 


IN 7 C 28 20 




CD TD it 2 


1 NT 02 8 40 


161 


IF ( RAI 3.NE.R 3IR i GO TO 1*2 


I N 7 0 2 Ii 5 0 




Rt$l(M=R3l8 


IN T 02 FcO 




PE 52 (K ) = 


1 N T U 2 o IJ 




1 F ( K. bw .2 ) o: TO In 


INT C2 t EC 




K-2 


I N YC2ci40 




GO TD loJ 


I NT J2 9CJ 


C CH EC * PC? 1 Mc^cCIpT £ RICFiT I NT ER SECT I "IN , EV M DATE IF REQUIRtD 


IN TC2S ID 


IE2 


IF ( 1 CF ♦ NE . IC) GO TJ lo2 


IMYC2S2J 




I S R= 1 


I Ml U24 -»J 




ou TC 163 


IN TC2S4C 


162 


IF 1 ( R A J R - R 3 J 3 ) /( P M P-ROIR ) * G T . 0*0) GC TO 163 


INT 02930 




AA=( AP<:FJ-A8 2R)/ (RAJR-RA K ) 


I » J T \jZ 960 




EA = AR2R-^A*rw I K 


IN T C 2 S 70 




AB=( AO *8 J-AT^R ) / (RB J 3 -kR IR ) 


IN To 29 80 




te = AR2 8-AB*F:MR 


INT C 2 S 0 




Rt S I A1 / ( A a— At) 


iNTGiOoJ 




P ES2 (K ) = A^RcrS I (K )+BA 


I NTOJO 10 




IF(K,tL*/i uC TO 155 


INT C 3 C 20 




K= 2 


INTOJj jJ 


C RE- 


INIT I AT E 


IN TOGO AO 


163 


a k2 R J = m F2 R 


INTC3C5C 




R 4 JL = p A I L 


INTODOtO 




RAJ R= P A I R 


I N 7 03 070 




R d JL =R G I L 


INT C 30 EG 




R3JR = P3 IR 


I N T 0 3 U vG 


140 


C .INT Ifji.l £ 


INT031U0 




UKl 7 E ( 6 ,210) 


IMTOJ 1 10 


210 


F J p ^ AT ( 1 L JQ° 1^0 T H R J JoH , N 1 PCINTS FOUND’) 


INT J 3i2o 




GG T C 1 :> 5 


INTCJIjU 


19 6 


IF(RCSl( I) .EC. ;ccc. U oO rc 191 


1.NTC214U 




IF (RES 1 ( c) • E 0 ♦ 4 w v. C • J ) G) TO 191 


: IN TO 3 1 50 




l F ( R E 52 ( 1) *l0.9JjJ.J ) oJ Ti 191 


IMTCJ UO 




lMRc$2( 2 ) .E-j.^OuJ.U) uJ T3 19 1 


INTO 31 70 




GO TC 153 


INTC318J 


C EVALLATUN CF I fsTERSECTI C f\S WhEN CORVES /RE ALM 1ST TANGENT 


INTO 3 ISO 


19 1 


YG >0.0 


INT 03200 




CALL FcNFTS IN2fNL,X 3, Y 3, Z3 , FT 1 f Y CO, ARRES ) 


IN TC 32 10 




Ki^ABM A R *< c 5 ( 2 ) - A 3 R E $ ( 1 ) ) fZ 


IN TC322 0 




CALL F t N P T 6 ( > , li,A4|T4,Z‘t,h'T^,TG3,AR Rt S ) 


I N 7032 30 




82- A E 5 ( jRRtSl2)-^nSLS(i))/2 


IN T C 3 2 SO 




RES 2 ( 1 )= ( ( - a **2)-(^**2) + ( 4RP2**2) i / (2. 0*4R?2) 


INT0325J 




RES2 (Z ) = F F S 2 U ) 


INTC326C 




RE Sl(l) =5 J*T { Ai3S ( ( Ri**2)-( RE 32 (1 )**2 ) ) ) 


INTO 32 70 




R E S 1 ( z ) = - 0 l S i ( 1) 


INT C3230 


155 


PtTURN 


IN TC 3290 




END 


INT G 3 3 CO 




SUBRCUT I Jt PEN 3 T^( 1 X , N Y . X , V , l , M T , YG , XR ) 


I N TO 3 3 I J 




DIREFS1CN Y (4 0) , <X(4J), 2(40,40 ),X(40 ),XR( lJ) 


IN TC 3 220 




I R- I 


IN T 03 3 30 


C SEARCH FCR J '1 F LOfccK Y LINE 


INT033hJ 




CO 1CI 1=1 , NY 


IN TC3350 




IF m 11 .Gc .YG) oJ T 3 102 


I N T 0 3 3 Cz 0 


101 


CO .'IT IMJE 


INT C 3 3 7 0 


1C2 


J= I - L 


IN TO 3 3 dO 


C CDAPSc SCAT F IP ZE<<D e\S3 


I N T 0 j 3 9 0 


C FEXT t LiNtS tX t C'JT t C i‘l FIRST St AR CF INLY 


IN TC3SOO 




YC = ( Yg- Y< j) 1 /{ Y ( JO 1 -Y ( J)l 


INTC34 10 




2X (11= U I 1. J+ I 1-2 I 1. J 1 1 *YC+Z( i,J 1 


I N TO 34 20 




IF (2X1 l.tC.nTj ol TO 107 


INTC343C 




-IT=2 


IN T 034 ^ J 


112 


CC 1 Cj I = I L', iX 


i ntgg^go; 




Z x I I 1 •=( Z ( I i J*- 1 1 - Z ( i ,j))*YC+Z(i,J) 


INT L J X c C'j 




ZA-l 21 I i J * 1 1 - Z ( I - 1»J 1 1 *( Yj-yIJI) /( Y( J* l J - Y( J) 1 +ZI I — 1 jJ) 


I NTGg‘1 7J 




IF (ZX( 1 ) .tt.hT) uf 'C 108 


INTcD-Ytc' 




IF ( ( ZXi I - Il-HT) / ( ZXI 1 l-i-.T) .L T.0.0) GC TO 1J ♦ 


INT D.H49J1 


10 3 


CENT 1 not 

/ 


IN T Jj 5 001 
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^ U> ONj»— 



. 



FILE: SWVOR 



FORTRAN PI 



GO TO IU 
C CAICULAT It. -i JF 



PFNET’AriONS 1 m 



£ I I HER 



TP l ANGLE 



INT C J S 1G 
INTO JSZO 



10 A 


X^=X (1-1 ) » (X ( i J -X ( l-l j )* (Y G-Y ( J n/ ( Y ( J*l J— Y ( J ) ) 


INTC35JC 




IF ( Z M . EJ.H N G J T ) ICS 


INT 03 5 AO 




IF ( ( ZX l I )-.iT )/ ( Z i-PT J .LT .0 .o ) GO TO ICS 


I NTC3330 




00 = A t£ ( hT-Z v ) mac S (ZX ( l-l )-H T » 


1NTC J5cG 




>$= >( I- 1 ) + ( XM-X( 1 -1) J *A > J( ZX( I-D-HT) /OD 


HjTO 3 5 70 




CO TC 11 J 


INTC35S0 


1 C 5 


>s=>p«-(x(ij-xMj*Mdi(ZP -nr ) / 1 a us < zm -ht j + a bs ( zx m- pt j j 


INTCJ5t»0 




lC TC 110 


I N T 03600 


C LCjI Cal FLNtT-uTlCNS ACCLV0LAT1CN 


IN TC26 10 


iu 7 


>A ( l ) = x ( 1 J 


INTO 3620 




1=1 


I TO ^6 30 




GJ TO 113 


INT u 2o AO 


10E 


xR ( P )=>(!) 


IMT 0 io 30 




GO TG ill 


IN T03660 


11C 


>R ( I R J = X S 


INTOic 00 


111 


IF < P .to. 2) GO TO IP 


I f J TO 36 oO 


113 


I K =2 


INTC36S0 




IB=P1 


IN fu 2 / OO 




GC TO L iz 


I N T 03 7 1 J 


116 


IF (IP. tC.U GO TC 117 


ANTC2720 




XR ( 2 J = 1 C G C . u 


IN TO 3 / 30 




GO TO I 14 


IN TO 3 / AO 


117 


>K C 1 ) =1 000.0 


IN TC3750 




AR ( 2) — l CCC. C 


INT 03 ?LJ 


114 


RETURN 


IN T03770 




ENL 


IN TL 3780 



SUBRCU INF l UPL7CX , Y ,Z ,X ,YG , f\X, NY, zPc S) 

C. LINEAR INT cP f CL« t IC'm CN CALIF- iT 1 IN S i J - F A C fc TO EVUJaTE 
C CP VALLE / 7 CO- KlM.Nm Tls XotYj. TnE RC SUIT PETURItG IS ZFES 
D I Rfc NS 1CN a (ru) , Y ( mu i ,L ( RO ,40 i 
C SEARCH PCh 1 C P LtfT X LINc 
JO 1 1 C = l » J X 
IF (X ( IC) .Ll.XGJ GO TO l 
CCMlNLc 
1= IC - 1 

SEARCH FCP J HP LOWER Y LINE 
DO 3 JC =1 , l> Y 
IF (Y (JC ) .GE . YGl GO TC A 
C C AT INJt 
J= JC- l 

AL=(Y(J + l)-l'(J)J/(X( I ♦ 1 J- a ( m 
EL = Y iJJ -AL* X ( l > 

7CR = ( iL*^) *• AL 
IF (YC.EO.Yi»R ) L) T) 7 
jJ T 3 5 

L'JwKtR TRIANGLE 



.R) 

I J 



IF ( V kj . \j 1 • Y 

C RESULTING L IS 
X3 =X (1*1) 

V 3 - V { J) 

Z3=Z ( 1+ 1, J J 
bt T C 6 

C RASllTINu l 15 INI 

5 X3 = X (1 ) 

Y3 = Y ( J+l ) 

ZJ=Z ( 1 1 J M ) 

C GENERAL CmULAT ICN: C 



UPrER TRIANGLE 



:CD FOR e-OT F TR IANGLES 



>l=X(Ii 
Y l = Y ( J ) 

Zl=Z(I , J) 

X2=x( i i ) 

Y2- Y ( J ^ I 1 
Z2=Z< 1 *L , Jti ) 

-AP = ( ( Z i-Z J) *( YI-Y2) -( Z1-Z2 )* ( YI-Yj ))/((< 1-X3 ) * ( Y l -Y2 ) - ( XI -X2 ) 
l*(Yi-Y3 J ) 

cP=<(Z1-Z2)-mF*( X l -X2 J J / (Y i-Y2 ) 

C P = Z I— A ? * X i- <P*Y 1 
ZR ES=AP*Xl vj+C^ 

GO TJ 100 

RESULTING Z IS 01 OWUING LINE 



I N TO 3 7 V J 
I N T C 3 E C C 
I i TO id 10 
I NTCJd2J 
IN TC i 6 20 
IMT0384J 
INTO Jd 30 
IN TC3fc60 
I NT 03670 
INT w3d 3 J 

nrcsijso 

I NTU3900 
INTC3SI0 
1 S P G 2 7 20 
I N T 9 3 J 
IN T L 3 5 AO 
IMT0395J 
I N 7 0 3 9 6 ) 
INT03S70 
I NT 0 39 SO 
INT C 25 SO 
INTCaulO 
I .NT J AO 1 J 
INTCAC2C 
I N T 0 A G 2 J 
I N T 0 A 0 AO 
INTCaCSO 
I NT 0 A 0 o J 
INTCAC7C 
INT04U9 0 
I NT J AO SO 
INTC41C0 
INT C A 1 10 
intoai^o 
INTCA13C 
IN i 0 A 1 AO 
IN TCALdJ 

in rcsico 

INT0A17J 
INTCaI EC 
INT0A190 
I NT OA^ 00 
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FILE: SrfVOR 



* 



FORTRAN P L 



N P 



S 



7 ZRES = (Z(l+l,JM)-ZU,J))*< Xu-X'm l/IXlIUl-Jdlllt/II.JI 

100 RtTUR'" 

E NO 

SLHR )CT I NL * 1 X A A Y (Z . ZMAX.NX ,NY 1 
C EVALUATION "’F 'AaHUM :P JF CALIBRATION SURFACE MATRIX 
DIMENSION il'tO.'.O) 

00 1 1= l, JX 

co i j = l ,-iy 

1 IF(Z< I , J) . ,t.zMAXJ ZfAX=Z(I t J) 

RETURN 

ENC 

iURRfUI I NE RI i\ *F V (Z , 2M1N , MX, NY ) 

C EVALUATi''N TF M'JMijm CP 3F CAL I RR A TI ON SURFACE MATRIX 

01 FfcFSi CN il-*J .*+0) 

00 J I- 1,0 a 

CO 1 J=i, IT 

1 1F( Z( I , J) .LC .Lin M ZMN=Z(I,J) 

RETURN 

£NC 

C END OF Sk>VJR 



INTCa 2 10 
INT0-.2ZJ 
IN TCAz AC 

I T T 0 A 2 A J 
IN T CAz AO 
IN I C A 2 tO 
INT J-,2 10 
IN JO ,23) 
IN i 0 ■» 2 0 
I NT 0 A i 00 
IN TCiilU 
I N Y C A A 2 0 
INTO'* j AO 
IN TC A3 AO 
INT 0 A a AO 
IN lOAioJ 
INT C ha 70 
INTJAJSJ 
INT Ci3SC 
INT Os AGO 
I NT0AA10 
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APPENDIX 3 



LISTING OF WV DR 



L?/Cfc/iS 
FILhJ i"VL n 






? 



S 



K^TMi\ FI U 



C 5 


TART OF »rVlN «'u** ‘ 




IN 7 000 10 




0 1 V E NS I N ilrllJli’Hl! lu) , * \ ( 4U t Y* (AC) , 




INT JJ02O 




lZA ( + •)»♦»)) v < r! ( 4 J I » V J l*i J ) t ZH (*tOi 'tc ) » X Aa( j ) i y \ \( *i 0 ) » - 




I J TClC iU 




2 L A X ( fXS/i^tC) ? Y i X ( H 0 i » Z 3 X ( ^ 0 f ^ i) ) 




INT J JO VO 


119 


l-CRMT (C j , , 




INT C t C 5 0 


131 


M'.^KF/.I) 




I N T Jjo oO 


1 JO 


F P 4 A T ( l j * L ) 




I NT 0 J 0 / J 


120 


( f-c .5 ) 




IN T C u u t o 


iJ 2 


FJ <«*Ol ( T 1C. 5) 




IN V C C 0 S J 


13i> 


F .) r * ( F I j . 5 ) 




IN T OuluO 


13 t 


PC ^M r (f"I/.3J 




INTtCliO 


C READ irt n ■ i ri UATA 




IN TOO 120 




Pc AC (2, 1 i\> ) IA, ;Y 




1 U TO 01 JO 




C'l L 21 i =1 , IN X 




I J 1 C 0 i i u 




Rt MJ( <r, 1 ;J ) aA( 1 1 




INTO J I jO 


1^1 


C C NT l N JL 




i.ncci lc 




0 1 12* 1 = 1 , si Y 




INTO 01 /O 




PeAL’(2,iJj) Y A ( l ) 




I N 7 J J 1 ■*» J 


122 


CCN T i N J c 




I 1 1 Y l I J 




l) j j 1 = l , N x 




IN' T 002 JO 




L -l 123 J = 1 , It 




iMuzn 




fl E A J i 2 , 1 2 0 1 <. A < j , j ) 




I h T 0 J ^ *_ J 


U 3 


cinr'Jct 




INTO 0 2 JO 




CJ 22 1 1 = 1 , lx 




I iN T Z J 2h- J 




RE.»J( 2 , 1 3u J >C ( i ) 




INTCt 2 ^U 


22 1 


C "l : 1 T If- JE 




l »N T j j Zoo 




DC 112 i = l , N Y 




IN T C C2 7 C 




RE AM ( < , liu) YC ( l J 




IN T J 0 2 0 J 


222 


CCM IN IE 




IN T 0 02 90 




CL) 223 i = 1 , IV a 




INTv.C3CU 




Of! 22 3 J = 1 , • J Y 




IN TOO j 10 




P E AC (2, I 2 j I Z3( 1, J I 




III T 0 0 J 2 ■ J 


22 3 


CC '1111,12 




In r j j o jo 


C TRANS! IV* ^LH’UIJN " 1 A 1 RI C2 3 




I f 1 1 J J J ^ 0 




CO 221 * = 1 ,l\f 




: INTCC35J 


22A 


X\/( l ) - Y U ( 1 ) 




INT CO: tO 




C ) 1 = 1 , U 




I N T C 0 j / vj 


225 


Y A x ( J ) = X A i i | 




I M T t C 3 c 0 




01 22 c I = 1 , J Y 




INT 0 0 3 V J 




DJ 2 26 J = t , 1 < 




I N T C J 4 00 


22t 


l\ >( I = 'M ( „ ,1 ) 




riTCtA I J 




DO 2 c 7 1 = 1,17 




I NT 00 ^2 ) 


227 


XbXiI)=YElii 




IN T C U 4 3 0 




JO c tb I = 1 , 'M X 




I N T uVn 0 


228 


YRX( 1 )=Xfc ( 1) 




I NTCJVNO 




CiJ 2 2A 1 - L , N Y 




IN T CC4F0 




L) J 2 * Ei J = 1 , 1 a 




I 1 i J 0 4 T 'j 


22 9 


2 u X ( 1 , J )-/(.), 1) 




I NT 0 ».) 4 o 0 


C FRCBt SclIi^S 




INTClASu 




YP A- C .0 




IN T J05>u0 




YPl<=5b. 




IN TO Ob 10 




XPA=C. 0 




IN TCG520 




XPO 2 C . 




I NT jo 5 JO 




XPD = -2 J . 




INTCC^AC 


C PR CG S A y Cu.'.S 1 '• JT 5 




IN TOO 0 30 




YrLC=-jO* 




I DTCJS jO 




Y R l P -30 . 




INTCu570 




YR I N = - ju. 




I>N T Job JO 




X P L C= -J 0 • 




irjTOub^j 




>R LP =3 0 . 




I N T 0 J c 0 0 




XR IM =- ^ 0 . 




I N 1 JvjoIj 




I C C P 3 - i 




IMTCL62C 




, h J 0 u P S = 2 




l N T C C c 3 0 




NO S l I 




INT0J6V0 




1 1 r -i 




IN T L Ct 3C 




I SCAN* I 




IN TOtcoO 




PLL.XP5- J. j 




INT j ot 70 


C EXPhPi HM 5 iUi;TiC\ 




I I TCct tC 


19j 


CM 1 iT'LE 




IN TO j6 90 


40 1 


f Cr-NAT ( r l<> j 




INTCC7C0 



12 



FILE 


: WVuR FORTRAN PI N 


p s 






RE AO 1 5,4011 PT 




n f 00 7 1J 




WR IT t (6, ~i u i 1 PT 




1NT Jj7dU 




If iNFS 1 w . C'j.2 I uC TO 500 




IN 7 C 0 7 30 




READ! 5, ^ C 1 i PS 




lMJo 7SJ 




WRIT E (oi'tvii i ° o 




IMCj750 




RcJulS.iCl 1 sLf A 




IN T C 0 7 c 0 




*R I T f ( _, <Ui ) ALFA 




1NT Ju 7 7 J 




FEAC(5.t01J P F 1 1 




mcc 7c j 




1*R I Tt ( b , AO 1 ) Ftll i 




i nt u/ so 




CALL i J TPL T ( aA, It , Z A , ALF A , Prt I I , 0 X , '1 Y , C 3 A ) 




1 M TJ JS0J 




A If C = A L F A-A P 0 




i .'ll L 0 d 1 0 




CALL I 0 f PL 1 ( aA , i A ,ZA ,ALFC , PH 1 1 , NX , NY ,CPC ) 




IN i 0 u ir 1 2 J 




alfc=alf;-xp- 




1 N TCoC jJ 




call iMf-L T (Xh , Y A ,za , alf J, Pul I ,\x , NY , CR3 ) 




i M 7 C t e S J 




PH IP=PH 1 i- <? 




1 N 7 J 0 b D J 




CALL I M T U L T UP, Y t .In, ALFA , =>HI3 ,Na, MY , C=> ^ 1 




IN i oUO o0 




w R I 1 E ( t> , '* Oil ) CPi ,CPc,CPL,oF0 




I N 1 > j d 7 0 




F UMf T< • T-" UL CPb »» tr lo. -j) 




1 NT OOOdO 




FJYI'.-FT-FS 




IN TCcciSJ 




PPA=C PA » F y Y N+ PS 




lNTuo^ 00 




fp^cpex 3 ;! no$ 




IN TO jn ij 




FFC = CFMfJi',\*P3 




IN f CCS3C 




PP0 = CP JXP^Y It-PS 




1m TO 0^30 




FRITtIo , a 1 J 1 FcA,FF£,FPC,RPD 




I 'j i L S 4 J 


403 


FORMAT l» HULL Pr S’, 4F10.2) 




I N T 0 u 9 


C RE f C if 'I') £ A 7 A 




I N TOO 060 




IF (NCSU'.L ..1 ; Tfj 501 




INTCCS70 


50C 


READ! 5,1 it J '■o.v 




1MT 0 JNdO 




WRITE (o, i jo I 33 A 




1 N TO OS V0 


* 


REZiO(5,135l RFo 




In rc loco 




WRlTEIo, iiol pph 




I NT J 10 10 




Rt tC (5,1 3 PI RFC 




1 J T C 1 C S 0 




WRITElo.lJc) Pr'i 




INTO lo iO 




REACI5, 1 j-j I ppj 




I N T 0 1 0 t J 




t* R 1 T E 1 6 , 1 i > P u 




IN To 10 tO 


iOi 


P.EaJI 5, 1 3 5 i Piil 




1NT0I060 




wRITE(o,IjPI Rjl'i 




INTOIO 70 




IF ( I 5C A i.L J.l l oj TO 181 




lf;T0 10-io 


18 C 


IF I I bF AN •O'I.ICjj) jU TO 30 5 




1 NT 0 10 SO 




1 1 T= 1 




IN TO 11 CO 




ICCPS=I 




INTCU1G 




15 CAT = ISCA m + 1 




INTO i 1 A) 




FSlA=F3iFr5. 0 




1 N T 0 1 1 J 0 




WRl TE ( t, lit! Pol. "I 




IN T 0 i 1 SO 


181 


CCM If' Jt 




IN TO 1 1 y j u 




F S- P S I X 




IN T 0 1 UC 




PaC=FS 




INTO! I 7 j 




FS T = PS 




1 N TO 1 i o 0 


C CALCLLXILS PRESSURE CCEf F IC 1 1 MS 




IN T J 1 1 bj 


30 0 


C0..T I' ' JE 




I r>j i ) i c. oo 




R R l T l < > ,N .j 1 llT 




IN TC 10 10 


SAC 


FORMA Tl » C* , 1 1 IC'JIIL'A Hj. ',151 




!NTv,i2du 




F0YN=PT"Pj 




1 TJ T J i d 00 




CFA=< FFA-i- i)/PjY.v 




IN TC UhU 




CP0=(Pp l-’SJ/Pji.. 




1 NT )Vl 50 




CFC= <P 3 C-PS J / ^ 




1 7 T 0 1 do J 




CPDM • VJ-FSJ / F . J Y . 




1 J f C 10 70 


906 


«R I TE l o , ,1J rr U » CP 3 » C PC » * °0 




I N T 0 1 x_ o J 


912 


F C k X i~ 1 ' oLCJ ol v S',AfiO. : i/ 




1 J TC 10S0 


C CORRECTS Pj -♦j^L'IPriwU T. £,.>0<-.c CP A- t IN C A1 I 6R AT l'l\ RANuE 


INYO i.JJ 




IF (NEC JCS.LN.nJ U T 3 jJi 




INT01 j 10 




ZAFAX--1C000. 




IN Tt 1 JAJ 




.CALL " \ X A . Y 1 ,Z 4 ' A a , IX ,' Y ) 




INTO 1 j 30 




IF ( C 3 A . C T .Z ' 1 a L | ’ON.Jl' 1 ') 




lNTCl-f^+j 




ificpc.jI.u'.o i-nnopf; 




IN ft I 3 tO 




IF l CP 0 . O i i •> A A 1 P ,0 = 1 . o 1 -P s 




I NT 0 I 3 0 J 




ZcPAX = -l lo j ) . 




IN T w 1 J 70 




CALL *U a \P 1 l Z-v , Z c UX ,\X , NY J 




IN Y 0 1 2 1 0 




IF(C°E.Ci.ZilAXj ■> )C : i . u i 


, 


I N T 1 1 o 




ZAM X = 1 0300. 




IN TC IS CO 
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FILE: UVCP 



fortran pi 



5 



*4 ? 



CALL MI 4ARY(ZA.ZAMN,f\X,AY) 

IF ( C PA .L f . Z A 1 1 4 ) PSC = u.99*P$ 

IF (CPC.LT.iAU4> FSC = 0.99*PS 
IFICFO.L l.ZAIlu) P SC =0.99* PS 
ZB^ 1*4- 10 jjj . 

CALL f 1 N A F y ( Z e ? Z E v i.M »NX , NY > 

I F( CP" .L r .2 3 'UN 1 P5C = c.99* PS 
IF (PSC . EC .?ST > o C T C 301 
ICCPi-lCCFSM 
96 3 wR ITE(o, >6 *) I Cl PS 

964 FCrNAT (• 1CCPS tS‘ , 15 > 

IF I lCCPS.uT. 1C) oC TC 180 

FS= PSC 
F S T = P $ C 

92C WRlTE(o,S2l> PSC 

921 FOPUTl* ^St IS ',F2J.5J 

GU 10 300 
301 CONTINUE 

C CALCLLATLS T WD ALP E PpI ANGLES FO° < I < £ 'll' 

CALL i N I cC S( YPa , riM , Y.RcJ ,YRUP, YRI N, 

1 C P A » C P ri , AL F , P F 1 , * A , Y A , Z A , 

2>B .YB.ZtU XA<, YAX ,ztX ,X cX ,Y EX ,Z SX » NX , NY ) 

ALF 1 = alF( l ) 

AL F2= AL F (2) 

PHI l=PdI 4 1 ) 

0 H 1 2= P H i ( i. ) 

C CALCULATES TWC ^ LF E Fni ANGLES FOP '1» 6 'III' 

CALL i 4 ISC8< XPA ,/?Z , >R LO »X PUP , X AI N » 

1CFA, CPC, Pr 1, 'LF.aAX ,YAX ,Z AX, 
2>AX«YAX,Z^X,XA,Y**iZA,XJ,YA,ZA,NY,MX) 

ALF3=ALf l 1 J 
AL F4 = A L f (2 l 
PHI 3 = PHi ( U 
Fh 14= PH I 12 J 

C SELECTS Trie PROPER ALF E Pul ANGLES OUT OF THE FUR VALUES 
FALF = ( ( ALFU- \lF2 + AlF3 + aLF4 J-AMAX1 ( ALF1 ,ALF2 ,ALF3, 

1 ALF4 J-A U !1( Al Fi , AuFc, alFj.A LF '.)) / 2 
FPHI M iPni l *■ P n 1 2 ♦ Pn i 3 +Pii I 4 ) -AM J,< 1 ( PHI 1, PHI 2, PHI 3, P F 4 4 ) 

1 — AM lNl(PFII,Prtii.|?UI i , PH 1 *») ) / 2 
900 V*« IT E (o , )0 1 ) FALF , r~>h 1 

9C1 F JRPA T( • SLLcCTEO' .2F10.2) 

C CALCULATES NEW Pj fFH DATA JF POSITION D AND CORRECTS PT 
ALFON=rALF-XPJ 

CALL I NT FIT ( XA, Y-> ,Za ,ALFu 4 ,F PH 1 , NX , NY ,C PO) 

ALFCN=F»LF-X°C 

CALL INTELT(XA f YA,ZA,ALFC J ,F PH I f NX , NY , C PC) 

909 WR ITF <o, 510) C^C ,CPO 

9L0 FCRMTf N E V C’C i C J C C AL LE 0 ' , «.F 1 5. 5 ) 

PSN= (PPC*C PJ-PPJ*CPC)/ (G FJ-C PCI 
907 wR ITF(;> , ,1 l ) fAu ,fPHI ,P3 . ,■> T 
9L1 F C E F A T ( • ITl-.JT iCN VAL'JcS' ,2F6.2,2FL0.2) 

C CON VEk GEM t <EilG 
950 CtLPS=PS4-PS 
551 V.KI Tt( j,952 J FS ,FS I ,CclPS 
952 FORM All* tPS^S V A LUE 5 ' ,3 F 15 . 5 J 
EP5PS= AS 5 ( ( P VI-* 3 S )/P S I 
So C FRllEIoi’iOll tPjFS 
961 FORMAT 4 * tPS^S I ' i r 2 5 . 1 5 ) 

EPSPSo=. JO JJoOOUoG99 * (PT**2 .G )- .UJ000098 *PT 
IF( ItPji.fcN.2) c P Sp SG = . L CG0CC0C009 *< PT**2. J)-. JO DO 003* PT 
IP(ItPji.E').2l t_» ^' 3 SL-=.CoJ\j0oJLGj9 ;!r (PT* t: ‘' s 2.jj'”.CGv.U\.CE* : PT 
IF(EPSP3.cT.EPSPSt» GC TO olO 
-GO TO iJi 
310 CONTINUE 

1 1 1 = 1 1 r + 1 

IF (1 IT.G'f. 1J) w.D TD 205 
FS=PS+RELX F 2* (PSN-PS > 

PSC=DS 
PST= PS 

524 RR1 It (6 » 525 > PS 



IN T U 14 10 
1 44 T 0 1 4 2 J 
IN 7c 1420 
IN r01't4o 
IN T 0145) 
IN TC 14C.C 
I NT 014 70 
I NTC I4EJ 
INTO 1490 
IN TO 15 00 
IN TC15 1G 
IN TO I d 20 
IN TO i. 5 30 
IN TC 1540 
INTO 15 50 
INTO I8 60 
IN T 0 1 5 70 
I NT0158J 
I N TC 1550 
INTO loOJ 
IN TO 16 10 
INTC1620 
I NT 01630 
1 N Tu 16 40 
INI Clfc 50 
INTO 1660 
INTC167C 
IN TO 16 00 
INTC1690 
INTO 1 / CO 
INTJ17 IJ 
IN T Cl 720 
INTO 1720 
I 44 TO 1740 
IN TO 1750 
IN TO 1 760 
I N TO l 770 
INTO 17 80 
INTO 1790 
INTClfiCO 
INTO 18 10 
I NTJ 1820 
INTCL830 
INTO 1840 
IN T01850 
IN TC IStO 
INT0I870 
INTC188J 
IN TO 1 850 
I NT 01 900 
IN TC 15 10 
INTO 1920 
INTO 19 30 
INTC154C 
INTO 19 30 
IN TC I960 
INTO 15 70 
INTO 19 80 
INTC15S0 
INV02UC0 
I N T u2 0 1 0 
IN TC2G20 
IN T 02030 
IN TC2C-tO 
IN TO 20 50 
I N T 02 060 
IN TC2C /C 
INT020 30 
I NT 02090 
INT02 100 
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> 



Flit 


: VVGF FORTRAN Pi N P S 




uy 

<N 

LT* 


FORMAT!' CU<i\ PS IS • ,F1 5 . 5 1 


in roil io 




GO TO POJ 


lN T G2i20 


303 


hrI T t (6 ,1 /l ) FaLF.FPH, PS'I, PT 


IN TC2 1 JO 


171 


FJRMA'I ( JF 1 a. 2 iF 1 5. 3) 


I JV02 140 




cj ro iv 


I N T 0 a. 1 30 


3 C 5 


Wft I 16(0, 30 7 ) 


IN TC2 RC 


JO 7 


FORMAT! 1 vCl'HEJ 1J0 TIMES') 


IN V 021 7J 




GO T 0 1VJ 


IM T J 2 i E 0 




END 


INTC2 ISO 




SU3R0UT I'Jc I iToCil iP. P1,ARP 2, a1L0»4RUP,4PIN, 


I N TO 220 J 




1HT1 ,HT2 . -<c3 1 , PcSa ,K*,y i, a! , 


IN TC2210 




2/2 ? Y l • 2. 2 »X2fYJfZ,JfX i t|Y4»Z4 * N 1 ? N2 ) 


IN T 0 2 2 ^0 




01 V EN$ i ~N : l 5 i ( 1 u ) , : £:>2( 10 ) , VRkE >( 10) , a1< 40) f Y 1( 4 0) , 


INTU22 JO 




1 ll 14 0 |4C) , X2 (4 J ) , Y2 <40 > , 22 A 4 ) ,4<J ) , X3 <49 1 ,Y ->i 40 ) t 


IN 7C2c4G 




2Z 3 ( 4 C » *C),A4<4U),YV(4C) *24(40, 40) 


INT 022 3 J 




RES1 (1 ) = >00 ) . 


I N T 0 22 oO 




PE Si (2) = 9 0 *j 9 • 


IN T 0 2 2 70 




$ E 5 2 ( 1 i = 0 j J o . 


I NT 022 JO 




Fc$2 (2) - 9 J 0 J . 


INTC22S0 




I SL^I 


I N V C 2 i 0 0 




ISR = I 


In IT 02 3 10 




0 A R = 1 . 


INTC2220 




AR 2R=AR ) N 


INT 02 jJO 




GO TC 13c 


1(4 T 0234 J 


15C 


*R2R-=AR2h + laR 


IN TC 2 3 50 




AR2PJ=AR << 


I N T 02 0 60 


152 


AR^l=Ar*2F—3r«Fl 


INTC237C 




AR^2=A^ 2R- ;kP2 


IN 


C CHECK FC'r- LC«K CALI FRaT 10 1 RANGE 


I NTC2 j90 




IF ( 7P"l . l T . 4PLQ) CO T C 150 


INTC24C0 




1F(APM2.LT. V.lJJ GJ TJ 150 


IN T 024 10 




K= I 


INT0242U 


C CALCIlATlS INITIAL P cNcT F« AT I CN PC1MS 


R ; T C2^iO 




call PtOPlM 11,0c ,a1 ,t 1, a ,rlTl , ARMl , ARSES) 


I N T u 2 h 4 0 




FA JL = A.<R£S 1 i ) 


* 1 N T C 2 4 5 0 




RA JK=APP6 S(2) 


1 N V C 2 4 c 0 




CALL PtNPTo ( 1 1*1 i,X2,Y2,Z2,HT2,ARM2,ARRES) 


I NT 024 (0 




RdOL«AFFLS (1 ) 


IN T C2 4 60 




RB JR= AS.P.c.3 ( 2) 


1NTJ2h )J 


C CALCULATES SUCcESIVE P EN ETRAT I fJN &C1NTS 


I N T J2 3 JO 




DO I A 0 i-1 ,10000 


IN T 02 5 10 




ARcR = AR2R»'JAR 


I NT 023 20 




ARM 1= a R2 -i- 1 


Rj T 02 5 30 




ARM2 = A‘<2R-4<=2 


IN T 0 t 5 40 


C CHECK FOR u:>>ER calibration ran^e 


I NTO 2 5 5 J 




I F ( A S'' 1 . uT . 3 5 IJO j GJ TC iV5 


I NT 025 >J 




IF(ARF2.jT.h"UP) uC TC Vjo 


I N r C 2 5 70 




CALL Pt'lPT j ( '1 1, l4,Xi, Yi,2 i , H T 1 , A R M 1 ,ARRES) 


I NT 02 3 RO 




R A I L-APPE3 ( i ) 


IN TC25S0 




RAI R = A,<ics( 2) 


IN V 02600 




CALL FL'IPfS ( 11,-Ja,X2,Y2,Z 2,HT2,AR 32,ARRE S) 


1 N T 02 <a 1 J 




FdlL=ARScS(l ) 


IN T 026 20 




RB I R= ARR E S ( 2) 


IN T02O30 




IF ( 1AJL . EM .10 JO . ) IS L= iO 


INT 02640 




IF (REJL.t ..l ) 0 0 .) l S l = 10 


INT0265O 




IKRBJR .l J i. IcJJ. ) I S?=10 


INT J2 6o'J 




IF ( f* 7 J ^ * c i » • 1)00* i 1 j iu 


INTC2L /O 




IFIPaIl.E,. IJOL. ) 1 U-1J 


INT 02 LEO 




IFIRPIL .ED. IjOJ. i I S«-= 10 


I NT 02 6 90 




iFiPPlh.EC.lJ >0. ) iSP= 10 


IN7C27CU 




IF (PAIR .E'J. 1000. J i j ', = 1 0 


INT 02 / 10 


C CHECK f-CP Cl<tCT HIT J| LEFT INTERSECT I JN, EVALUATE IF REQUIRE D 


I NT 02 720 




IF 1 I SL. NE. 10) JO TO loA 


INT C 2 /JO 




- 1 SL= I 


INT ) 2 / 4 J 




CC TC 1 3 J _ 


INT C2 7 50 


164 


IF I « A 1 l. N t . R 31 L ) GO TC 180 


INT *J 2 tt J 




RES1 (n l^'UL 


I NT 02 / 7 J 




PES2 (K) = Ap.^R 


I M T C 2 / 2 0 




I F ( r<, . L. . 2 ) o- TJ 155 


IN f C2 7 90 




K=2 


IMTJ28 JO 



'b 
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FORTRAN Pi 

1 d 3 

i at RMEOI Af E l eft 



FILE: WVO° 

GJ TO 
C CHECK. ITk 

18 0 IF l ISL.Nt. 10) GO TJ 1-60 

ISL= I 
CO TO idi 

1 1 C IF (( F^OL^t JL) /( F A AL-R6I L) .GT .0.0) 
AA=< AF 2R J- n 2^ ) / ( K; jL-r.Al L ) 
EA=APcR-/A*7 AIL 
A3-( AK2 :< J-^2 k) /(FdJL-PdlL) 
BB^AR^-AN^FRiL 
KESl ( K) = ( eb-f A) / |/J./ 6 ) 

RE J 2 l k ) = A A * < L j l ( n ) + 3 A 
IF (K .LC. fc 2) G ) n i?b 



«m*i' * 

INTERSECTION , EVALUATE 



GC ti 16 3 



IF REQL1REJ 



183 



K = ? 
CHECK rJR 



01REC1 HIT 



IF ( l^F .M.1J ) GO 
1 5*M 



J i 
TO 



RIGHT 

iol 



InTEKSECTI CN.EVALUAT t IF REQUIRED 



GO TO 163 

161 IF <Rnrc.Kfc.9Bl R> CO T C 182 

KEil( K ) - r i i * 

PfcS2 (K ) = A~;2? 

IFlK.fcO.2J oJ TJ 155 

K- 2 

GO TO 163 

C CriECK FOR 1 inTER^EJI ATE F I oHT I NT ERSECT ION » EV ALJ AT E IF R E CU I R ED 
182 1FIISF.ME.1J) GO TO 162 

I S R-l 
GO TO 1 c 2 

162 1FUFAJF-^Ja)/IRMR-R6IR ) .GT.C.O) 

A A - ( £P2Ro-Ar\c P ) / ( FAJn-hAIR y 
8A= AR2R- AIR 
*6= ( AF2P J / IP3J^-RB n ) 

8d = AR2r<-A3*^IR 
RES 1 ( K ) = ( A A- A A J / ( A A- AQ ) 



GC TO 163 



b(2J- - <t SI l) ) / 2 

( -<4.***: )+l A.<p 2 * ; 



c) ) / ( 2. 0*A*P 



R2=A 3 S(AR‘ : 

FES 2 (1 ) = ( ( R 1 * *2 i- 
RE 52 ( 2 ) = R L 'iZ 1 1 ) 

RES 1 m = SO\T< A*n l *1**2)- I Kfc S2( 1)**2 ) ) ) 

FfcSi ( 2 ) i m 

155 -RETLkN 
ENL 

SUOR CUT 1 K h FLNPT S ( 4X , NY ,X ,Y , Z , FT , V G» X P ) 
CI^ENS I JiM Y ( AC) #«_>( ^ C ) U ( <iyj t 90 ) ♦ K( AG) , *R ( 1 0) 
I R - 1 

C SEARCH FOR J L F i C V\ E r V LINE 
un iol 1 - 1 1 r j y 
IF 4 Y ( I ) .Lc.Yu) jC TC lu *: 



INTO^J n 
I N T H20 
INTO 28 JO 
IN T J 2 6 *i J 
lNT02o:>0 
IN T 6 2 8 f; J 
I NT U2 b /O 
IN T C 2 ci c 0 
INI 02E50 
1 N T 025 00 
IN f C <_ S 1 0 
IN T 02520 
I 0102936 
IN T C 2 S aO 
1NT0295J 
INT C2S60 
lN T J2^7u 
I NT 02 5 3 J 
IN TC2 5S0 
I NT G 3 G ^ 0 
INTQJJlJ 
INT03G.r0 
INTuiO’iO 
1 N TO J 05 J 
INTO 3 0 50 
I NT 03060 
INTL3C/C 
IN T 03 0 80 
1NTCJ 090 
INTC31CG 
I N T 0 3 i 1J 
1NTC312J 
INTO 3 130 
INTG31 50 





RtSZ (K ) = AA*ruS 1 ( K ) *SF 


IN TC 3 150 




IF 1 K.t J. 2) j J TJ lD5 


INTJ316J 




K= ^ 


INTO 31 io 


C re- 


1 M Tl Me 


I N I C 3 1 c C 


163 




IN T 03 1 90 




«AJl = P^ 1L 


I NT 0^2 C J 




R A wK =-< A IX 


IN T C 32 10 




R B J L = P i 1 L 


1 N TO J 2 20 




i? 


i m T C J 2 3 0 


14C 


CONll.'JLc 


i* 1 7 0 3 2 90 




WR ITtlo , )) 


1 N TGa2 5 ) 


210 


FCF^ATC L'.'JP 1^0 7 J-hGU Jn, *n RIIMTS FTJMC'i 


IN T C 3 2 1 0 




60 TO 155 


I N T j 3 2 7 J 


195 


If (R£Si 111 .EO .9 JUJ.'J 1 JO T 0 19 1 


IUT03230 




IF 1 RcSl (2) .c J.9J00.0 J J 1 T ^ 19 1 


IN T 0 j2S0 




IFIRFS.m .l ). y^J. J) TJ 191 


INTO 33 00 




IF ( » fcS_ (2 ) . t J.906 J .0 1 o J ri 191 


INTC J3 10 




GO TC 155 


INTO 3320 


C EVALUATION )F riTER^fcCf I--NS when clrves are alnjst tangent 


INTO 3330 


191 


Yb C=0 . j 


M TC 33^0 




CAlL PC OPT J< -I2.WI ,X3 ,YT ,ZJ ,HT1 .YJOtAPPESl 


IN TO 3 3 50 




RWtSiafi:-. lSU)-a*;c:>< i/^ 


IN 703360 




CALL t- t.NF I ") ( \J * i j 1 , <-t ,Yh , l* ,nT^ ,Yu J, AF PCS J 


IN TC J 3 70 



1 NT0333J 
INT0J3^0 
INTO J4 00 
1NT0J510 
IN T C 3 5 20 
1NT CjA 30 
I N J J j550 
IN TC3550 
I N I 0 j 5 c J 
IN TO 3 4 70 
IN7C34£o 
1NT 0 j VJJ 
INTC35CC 
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F J IE: WVUfi 



FUR TRAN ?L 



IC1 


CuNTIhLE 






IN V 0 3b 10 


102 


J= i-i 






I NTU-520 


C CGAivSL SCAN fCh ZERC FYSS 




* 


IN 702*20 


C NEXT 4 LINES e EXtGLi'Ci) IN FIRST SEARCH CNIY 




i N T ) 3 5 '♦ J 




YC = ( Yt-r (J J 1 / l r < J ♦! )-» ( j )) 






IN Tv- 3330 




2X1 1 ) = Li ( l , J <•!)-.. <1 i.l i $Y C ♦Z ( I » J J 






InTU35uO 




1 F ( Z XI L ) .EQ.HT) uJ TG i J / 






I NT 03 3 70 




I li = Z 






IN TC 2 5 80 


112 


JJ 1 C3 I=1'5 ,I\X 






I NT C3 5 30 




ZX ( i J=( Z( 1 » J + 1 ) - Z l 1 > J ) ) * YC ♦Z ( 1 iJl 






iNTGjoGJ 




ZM=( i ( 1 , .u ) -Zl X-i, J J I* (YJ-Y ( J ) J/( Y( J U)-Y< J J ) +Z ( 1-1 , J) 


INTC3EIC 




IF I Z X 1 1 ) . n J. rt I ) U C rc I Gd 






IN 1 J j t 2 0 




IF I ( ZX( I- U -HTJ/ (Z<( I J-FTJ .LT .L.O) 


GC TO 104 


I N 1 6 3 3 0 


103 


C JN T I Me 






IN TC JC^G 




GO TO lid 






INTO 56 30 


C CALCULATION IF 0 EN ETF '7 1 DNS IN EITHER T R I ANGLE 




INTU3OO0 


1C< 


XM = > { 1-1 ) «• 1X1 i J-X ( i-I J )* IY u-Y l J j ) / ( Y U + iJ- 


Y ( J 1 J 


I M 1 C 3 1 i J 




IF(Z’l.cO.hf) GJ T 3 Ui 






I \ i OjOuJ 




IF I ( /X ( i )-FT )/ U' -FT l ,L T . J . J) G3 TO 


IC3 




I jN T C 3 6 3 0 




00 -A ►'•SI UT-i. v )+AbS (Zx |i -1 )-nT J 






INTbi 7 Cu 




XS-X ( I-1JMa<-aI 1-1) )^oo(ZA{I-i)-HT)/OD 




1NTC37 iJ 




G C T t I I 0 






I J 7CJ?^C 


1C5 


> 4+ ( *( ii -AM)* A3^( L» -HT ) /( AcM l'A 


— ri T ) AB 


> ( ZX ( I ) -nT J i 


I 'JTO 3 7 30 




GO TC llu 






IM 763740 


C LOGICAL Fti\ 1 1 5 A T IDAS ACCLFULAT 1CN 






IN TC 2 7 20 


10 i 


XR l 1 )= a( J J 






INT 03 7cJ 




1 = 1 






INTCd 776 




GO 10 113 






INTO 3 7 HO 


108 


XR ( IP J= X ( I J 






I NT 03 7 40 




bu TC 1 L 1 






IN T 6 38 CO 


lit 


*Rl I R)= >5 






I N T 0 3 o 10 


11 1 


IF C IP . tO ,2 J CO TO 114 






IN T6 38 20 


113 


I R=2 






IN T03d2J 




1R~ Ml 






INT 03 c 40 




GO TC 112 






: INT0383J 


lit 


IF ( 1 R.EU.l ) oC TC 117 






INT C3E tC 




XR ( 2 ) = Uoj • 0 






I NT 0 3b 70 




CO TC 114 






INTC3t EC 


111 


>R i I)=1COO.O 






INT U 38 00 




XR ( 2 )= 1 C 0 J ♦ J 


1 




I NT03900 


114 


RETURN 






INTC3S10 




ENO 






IN YC3420 




SUER7UT 1ME 1 ITPlTIX.Y.Z.xG, yG.NX, ‘Y,Z"S) 




I T 0 3 v 30 


C LINEAR IMtRFCLATJCN Cl C ALi AT i CN SURFACE TR 


EVAL JU E 


INTC3S4C 


C CP 


VALLE A i C j-O.xl) i ,N . Tc S Xj.Yj. Tn t RESULT RETUKNCJ IS ZFES 


INI 0 5^ 30 




GI W ENS i IN X (-.0) ,Y ItJj.ZUJi^OI 






INT03SL0 


C SEARCH FCk I LP LEr T X LINE 






IMTC5S 70 




00 1 I C = 1 » N X 






I Nl 0 3 C J 8 J 




1 F ( X i IC J .CL.Xb) oC TC 2 






I iM T C 3 S c . C 


1 


CONTI iJ L E 






IN T 040 JO 


2 


i= iC-1 






I NTGhO 10 


C SEARCH PCR J Cr LOWER Y LINE 






IN TL4C20 




00 3 J3=l, IY 






IN TO 40 30 




IF (Y (JC J . L E . Y G ) CO TO 4 






INT04O40 


3 


CCMIMc 






IN TC4030 


4 


J= JC- I 






INT 04060 




AL= (Y(J +1J-Y (J ) )/ (X l I + 1 ) - X ( I)) 






I N T 6 4 0 7 0 




BL= Y l J) -A L 5 * X 1 1 ) 






IN T C h C t O 




YC R= (AL *XG i + TL 






I N ♦ O^L 00 




IFiYG.cC.YLR ) uJ TT J 






I NTC41 CO 




IF ( YG.iT. YC U o J TO 5 






in ro^ i io 


C RESULTING Z IS I 1 LJ-^tR TR1AJGLE 






1 NIC f 120 




X3=X U + 1J 






m 6 4 I 30 




-Y3= Y ( J ) 






I.NT04 i-tu 




2a= Z I1U.JI 






IN TO Hi 50 




GO TO 6 






INTO 4 1 u 0 


C RASULTnc L IS IN UPPER TRIANGLE 






I NY 041 70 


5 


X3=X ( i ) 






I M TO H 1 80 




Y3=3( J+ 1) 






INT 64 1 SO 




2j=ZU 9 J'* l) 




, * 


INT J 42 OJ 
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FILE: W VGP F i)9 TRAN P I 



N P 



S 



c general calcilati jn :gc jU fgr dCTh Bungles 

6 xi = x ( i j 

Y 1 - Y ( J J 
Zi=Z< I 9 J ) 

X2 = x uu; 

>2 = >( J* L J 
Z2-Z ( hi, Jill 

AP= ( 4*1- <3 J * (Y l -Y ^ j- (Zi-ZZ )M YI-Y3 ))/((Xl-Xi)*(YI-Y2)-(Xl-X2) 
1*( Y 1 — Y -5 J ) 

6P= ( ( Li-U J - *P* ( X 1- XZ ) }/( Y 1- Yc ) 

CP=Zl-A P* <1 -~P*Y i 
ZR ES =A P * /G+HP* ru+C 3 
GO T n UJ 

C RESILTING L IS Ci\ OIVUII^ LINE 

7 ZRF3=(Z { 1 f I, j* i)-Z( 1 , J ) J *4 Xo 

100 PETUhh 

EN J 

SUM 0 CUT INE MAX ARY (Z f ZMAx f IX 9 
C EVALIATICh UF N/Xivj/' ^P CF ' 

DIMENSION Z(AG,Av.J 
CO 1 1 = 1, AX 
D1 i J-l , J Y 

1 IFtli 1, J I .Jt.ZMAA) ZMAX=Z(i,Ji 



X(i))/lX(ItI)-X(nifZ(l,j) 



ir ) 

CALI [JR AT 1 ON 



SO- FACE MAT R IX 



RETURN 

End 

SURRCUT I 'IE M [:UR Y (Z , ZM IN f N X,N Y ) 

C E VALLA T 1 CN Lr i v I A i V: J ^ CF CF CAL i 0 < AT ION SO 0 F AC E MATRIX 
DIMENSION Z I A C , A l ) 

CO 1 I - 1 , MX 
CO 1 J=1,NY 

1 IFIZ( I, JI.LE.ZMIN) ZM I N=Z ( I , J) 

RETURN 

END 

C ENC OF WVD* 



IN T C a 2 IJ 
I N TuA 2 20 
IN TC <,2 30 

i r j T o a l t o 

1 N T 0 >2 b J 
I N T L A * u U 
i NT ) 't2 / J 
PI7Ct2cJ 
1MVCA2SM 
IN 1 C a 3 o J 
I N Y L A J 10 
IN T DA j 20 
IN TO'tSiO 
l n Y 0 A 3 A J 
IflTiJ-ObJ 
I N Y c^3^ w 
IN T J A 3 
I N T C A A ci j 
IN T C <3 SO 
1 1* ‘I r 0 A A j J 
1 N T C A A i 0 
IMT Jv, 2u 
1 NT 0 A A3 J 
INTGAAAU 
INTCAAbO 
1 N T C A A oj 
IN T C a A ? C 
IN TO aA d J 
1 N T CA4 AO 
IN T CASCO 
INT’J'O l-J 
I N 7 U*t 5 20 
l N “f 0 *i b 3 J 
I NTOAbA J 
I N T C A 5 5 0 
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